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The Timoshenko effects play an important role in the 
analysis of structural components frequently appearing in 
the design problems* Keeping this in view, economically 
efficient analysis techniques are essential. With this 
aim, the elemental matrices for the finite element formula- 
tion of various structural components, viz., beams, rings, 
sectors and discs, have been deduced logically from the 
governing differential equations using the Galerkin method. 


The Timoshenko- beam has been analysed using two 
approaches. The first approach applies to a uniform beam 
characterized by a single governing differential equation. 
The second approach uses two simultaneous governing differen 
tial equations and is applicable to any beam. Though the 
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elemental matrices for second approach are given for a 
linear variation of the cross-sectional parameters, the 
evaluation of matrices for any other variation is straight- 
forward. 

The equations of motion of a ring or sector, with 
and without Timoshenko effects, are improved to include the 
rotational speed effects. The finite element equations are 
systematically derived from these governing differential 
equations. A quintic polynomial, satisfying the compati- 
bility of derivatives upto second order, has been used for 
the ring finite element. It is shown logically that the 
order of the interpolating function remains same, in case 
of the analysis of rings and sectors, with and without 
Timoshenko effects. 

In a similar logical fashion the finite element 
equations have been derived for a uniform disc, with and 
without Timoshenko effects. The results for a tapered disc 
are obtained by idealizing it with the constant thickness 
elements. 

It is very easy and straight-forward to satisfy all 
the natural boundary conditions in the analysis of structural 
components, when the finite element equations are deduced 
using the Galerkin method. This is illustrated by the results 
for the various problems discussed in the thesis. 
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The critical ratio of the inner to outer diameter is 
evaluated for a centrally clamped disc. The lowest frequency 
for ratios above'; the critical value occurs for the one 
nodal diameter and zero nodal circle mode instead of the 
zero nodal diameter and zero nodal circle mode. The change 
of critical frequency, with rotational speed, has been 
studied for a centrally clamped disc. 

It is illustrated that the order of interpolating 
function does not change when the Timoshenko effects are 
included in the analysis of structural components. The 
rotational stresses pose no additional complexity and cause 
no change in the interpolating function when considered in 
the analysis of rings and discs. 



CHAPTER I 


INTRODUCTION 

The advancement of technology and high speed machines 
and vehicles has necessitated the need for a more accurate 
analysis of the structural and machine components. The 
classical solutions of vibrational problem of various 
structural components, like beams, sectors, rings, discs 
etc., have long been obtained. These are adequately 
applicable to lower modes of vibration. The higher speeds, 
generally, result in excitation of vibration modes of 
higher order where the Timoshenko effects, namely, the 
rotatory inertia and the shear effect, play an important 
role. The governing equations, with the Timoshenko effects 
included, become complicated and the analytical solutions 
become difficult or impossible. Even for the simplest 
component, the Timoshenko beam, the analytical solutions 
are not possible for all sorts of boundary conditions 
and for all variations of thickness of beam. Simple 
geometric shapes and boundary conditions can only be analysed 
analytically. For other components it becomes virtually 
impossible to obtain analytical solutions. Keeping 
those problems in view, various attempts have been 
made to obtain numerical solutions to various structural 
problems. In the present study the finite element method 
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has been used to solve various vibrational problems of the 
structural components. 

1.1 LITERATURE SURVEY 

The most important and basic structural component is 
the beam. The error introduced by neglecting the Timoshenko 
effects is small for lower modes of thin and long beams, but 
turns out to be significant for higher modes of all beam. 
Archer [l] has given a finite element formulation based on 
classical theory. For a general beam, it is essential to 
develop a convenient analysis technique based on Timoshenko 
beam theory [2] i.e., including the effects of rotatory 
inertia and the shear. 

Huang [3] has derived the exact frequency equations, 
for simple end conditions, of uniform beams based on 
Timoshenko beam theory. The equations are not conveniently 
usable for design purposes and cannot take-care-of the 
variations in the cross-section. Moreover, it is quite 
cumbersome to derive the frequency equations for more 
involved boundary conditions and for beams on multiple- 
supports. Thus, the requirement for a numerical technique 
to analyse the Timoshenko beam. 

A number of Timoshenko beam finite elements have 
been proposed in literature [4-15]. Thomas et al [4] have 
very comprehensively discussed the earlier elements and 
have proposed a new element with six nodal parameters. The 
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elements have been classified into two groups, simple 
elements and complex elements. An element with two nodes 
and four nodal parameters is termed as simple element and 
the one with more nodes and/or nodel parameters is classi- 
fied as complex. Table 1,1 gives the summary of different 
elements and is reproduced from Thomas et al. [4], None of 
these shape functions have been deduced logically, instead 
the interpolation functions for the transverse displacement 
and the shear have been selected arbitrarily. In some 
cases, it becomes difficult to satisfy the natural boundary 
conditions. A detailed review of these functions is given 
in Chapter III. 

Table 1.1 

Summary of different Timoshenko beam elements [4] 


Element Reference 

No, of 
Nodes 

No. of 
parameters 

Parameters 

A 

[5-11] 

2 

4 


W, 

B 

[12,13] 

2 

4 


w, ^ 





dWj^ 

dWg 

C 

[14] 

2 

8 

®b’ar'" 


D 

[13] 

3 

7 

dw 

at each endj 





at middle node 

E 

[15] 

4 

8 



F 

[4] 

2 

6 

W, ^,Y 
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The sector or the ring is equally important structural 
component. The understanding of vibrational characteristics 
of complicated structures, such as gears, aircraft fuselags, 
missiles, etc., is facilitated by understanding the dynamic 
characteristics of rings on radial supports. Dynamic 
characteristics of free-thin rings have long been investi- 
gated by several authors [ 16-19]. Different methods have 
been proposed by several researchers [20-34] for the vibra- 
tional analysis of thin rings with and without radial supports* 
Rao and Sundararajan [22] have determined natural frequencies 
and mode shapes of free and radially supported thin and thick 
rings by constructing the frequency determinant using the 
differential equation approach. The method proposed is not 
suitable for large number of supports because the order of 
the frequency determinant increases six-folds with increase 
in number of supports. Moreover, even for small number 
of supports, their results are inaccurate [23,24]. For a 
thin ring with large number of equi-spaced radial supports, 
Sahay and Sundararajan [25] presented another method, which 
is not suitable for small number of supports. Also, the 
violation of inextensibility condition of the centre-line, 
assumed in their formulation, results in erroneous values [24]. 
McDaniel [26] used the complementary transfer matrix method 
for dynamic analysis of the thin ring on elastic supports, 
Murthy and Nig am [23] used the transfer matrix approach to 
investigate the problem of the thin ring on equi-spaced 
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radial supports. It appears that the adhoc assumption, 
that the vibration mode of the complete ring can be 
obtained by repeating the mode of one bay, has resulted in 
missing certain frequencies for the ring with more than 
two supports, Mallik and Mead [24] applied a wave approach 
to solve the problem of the thin ring on equi-spaced, 
radial supports* The pre-condition of the method is that 
the structure must be rep£titivo. This leaves the method 
to be of limited use. 

In most of the earlier works, the shape functions for 
rings have been deduced from those for cylindrical shells 
[27-29] by specializing them to two-dimensions. Sabir et al, 
[30-33] deduced the shape functions for the analysis of 
static arches from the shape functions for cylindrical 
shells [27-29] and made a comparison of these, Sabir and 
Ashwell [34] extended their work to dynamic analysis of 
thin arches/sectors and thin rings. The shape functions 
suggested are not general enough to take-care-of all types 
of geometrically different rings and sectors. Also the 
number of elements required for convergence of results is 
considerably large, varying from 15 to 20, 

The classical solution, given by Hoppe [l6], neglecting 
the effects of shear deformation and rotatory inertia, is 
applicable to low nodal diameter modes of thin rings. Care- 
fully carried-out experimental studies [l8,35] have revealed 
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that, even for thin rings, considerable error is introduced 
when the number of nodal diameters become large. Kirkhope 
[36] has reported that Federhofer [37] has derived an 
exact solution in terms of Bessel’s functions for free 
rings including the effects of shear deformation, rotatory 
inertia and extension. The results are not readily applica- 
ble to design problems [36], Seidel and Erdelyi [38] have 
used an energy approach to derive an approximate solution, 
with assumptions regarding nonlinear normal strain through 
the cross-section, and an average shear angle. The cubic 
frequency equation obtained is further reduced to quadratic 
expression by showing the effects of mid-surface extension 
to be negligible, Rao and Sundararajan [22] solved the 
differential equation of inextensible motion for the natural 
frequencies and mode shapes of the thick rings. The assump- 
tion of inextensibility results in a quadratic frequency 
equation. For large T/R ratio, the results of references 
[22,38] deviate considerably from experimental values. 

Using the energy approach and making certain simplifying 
assumptions, Kirkhope [36] derived a frequency expression 
for in^plane vibrations of a free in-extensible thick ring 


r ( 1+n^ Si— 2^+^ 1) 

i\ aAGR^ 


aAGR 


( 1 . 1 ) 
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In the design of turbines, nuclear reactors, machine 
structures, etc., discs are an important component. The 
analytical solution for vibrations of an unconstrained, 
uniform, thin disc was carried-out by Kirchhoff [39]. The 
work was further extended to include the effect of rotational iri- 
piane stresses in free and centrally clamped discs by 
Lamb and Southwell [40] and Southewell [4l], Mote [42] 
has proposed an approximate solution procedure, based on 
Rayleigh-Ritz technique, for the analysis of initially 
stressed circular discs. The analytical solutions are 
almost impossible to obtain for discs of variable thickness 
and for discs including Timoshenko effects. For such pro- 
blems, numerical techniques, like finite element method, 
have to be used. 

The dynamic behaviour of the discs can be predicted 
by annular and sector elements. For symmetric vibrations 
of discs, the annular elements are more convenient to use. 

The variation of displacement in the angular direction is 
exactly represented by trigonometric functions, therefore, 
the annular elements can be termed as semi-analytical. 

Polynomial approximation is used in the radial direction. 

Thus, the problem reduces from two-dimensions to one- 
dimension with a considerable saving in nodal parameters and 
increase in accuracy. 
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Several annular finite element formulations based on 
thin plate theory have been reported in literature [43-49] • 
Kirkhope and Vtfilson [43] have proposed an annular finite 
element using different cosine functions for each vibration 
mode to represent displacement function in angular direction. 
A cubic polynomial is used for representing displacement 
function in' radial direction. The element has 2 nodes, 

4 nodal parameters and linear or parabolic variation in 
thickness. The element was further developed by Kirkhope 
and Wilson [44] to include the effects of rotation and pre- 
stress. Kennedy and Gorman [45] used the original element 
to study the effects of centrifugal and thermal stresses on 
dynamic characteristics of variable thickness discs. There 
are many other formulations proposed for thin discs [46-49] 
in literature, 

A number of elements have been proposed for dynamic 
analysis of thick discs [49-53], Wilson and Kirkhope [50] 
have developed a 2 node, 8 nodal parameters element with 
linear variation of thickness. The nodal parameters being 
the displacement, radial slope, and the radial and tangential 
shear rotations, A 3 node, 6 nodal parameters and parabolic 
profile annular element has been developed by Hinton [51] for 
axisymmetric vibrations of thick arbitrary plates. The 
nodal parameters being the displacement and slope at each 
node. Mota Soares et al [49] have developed a 3 node, 

9 nodal parameters annular thick element with parabolic 
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profile. The parameters are the displacement, and the 
radial and tangential slopes at each node. An annular 
element for the axisymmetric vibrations of circular uniform 
plates has been proposed by Raju and Venkateswara Rao [52j, 
The element has 2 nodes, 8 nodal parameters, the displace- 
ment and its derivative, and the shear rotation and its 
derivative are the parameters at each node. The element 
is also used for nonlinear dynamic analysis of circular 
plates, Mota Soares and Petyt [53] have developed an 
annular finite element, for dynamic analysis of arbitrary 
discs, based on Mindlin plate theory. The element has 
2 nodes and 12 nodal parameters. The parameters per node 
are the displacement, radial and angular slopes, and the 
radial derivatives of displacement, and the radial and tan- 
gential slopes, 

A summary of various finite elements for analysis of 
thick discs/plates is given in Table 1.2, 

1,2 PRESENT STUDY 

■ It was observed through the literature review that 
even though considerable amount of work has been reported on 
the finite element formulations of various structural com- 
ponents, most of them appear to be lacking the logic of the 
finite element method. Almost in all cases the interpolating 
functions are selected without considering the compatibility 
and completeness requirements. The functions selected are. 
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Table 1.2 

Summary of various thick disc/plate elements 


Author( s ) 

Reference 

No. of 
nodes 

No. of 
para- 
meters 

Parameters at 
each node 

Mota Soares, 

Petyt and Salama 

[49] 

3 

9 

displacement, and 
the radial and 
tangential slopes 

Wilson and 
Kirkhope 

[50] 

2 

8 

displacement , radial 
slope, radial and 
tangential shear 
rotations 

Hinton 

[51] 

3 

6 

displacement and 
slope 

Raju and 
Venkateswara Rao 

[52] 

2 

8 

displacement and 
its derivative, the 
shear rotation and 
its derivative 

Mota Soares and 
Petyt 

[53] 

2 

12 

displacement, radial 
and angular slopes, 
and their deriva- 
tives 


in general, of unnecessarily higher order, thereby increasing 
the number of variables. Keeping in view, the importance of 
the analysis of these components and frequency of their 
occurrence in the design problems, it is desirable that the 
analysis techniques used must be economically efficient. With 
this in mind, the elemental matrices for finite element formu- 
lation of various structural components have been deduced 
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logically from the governing differential equations using 
the Galerkin method. 

In case of the Timoshenko beam formulation two 
approaches are used. Approach A, using a single governing 
differential equation for a uniform beam. Approach B, using 
two simultaneous governing differential equations for a 
general beam. Though the elemental matrices for Approach B 
are given for a linear variation of cross-sectional para- 
meters, the evaluation of the matrices for any variation 
of taper is straight-forward. Both of the new elements 
fall in the category of simple elements as classified by 
Thomas et al [4], The derivations are given in Chapter III. 

•i 

The derivation of the elemental equations for the 
rings and sectors, with and without Timoshenko effects, is 
given in Chapter IV, It is shown logically that the order 
of the interpolating function is same in both the cases. 

The interpolating function for a line element using compa- 
tibility upto second order derivatives has been used for 
the analysis of rings and sectors. 

Chapter V deals with the derivation of finite element 
equations for the discs, with and without Timoshenko 
effects. The derivation is given for uniform discs only, 
because the results for rotational stresses for a tapered 
disc are not available in closed form. Hence, the 
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matrices have to be obtained by numerical integration in 
case of a tapered disc. 

The efficiency and accuracy of the formulations have 
been illustrated by comparing the results for various 
problems with the available results wherever possible. 
Results for many of the unsolved problems are also presented. 
The ease with which the natural boundary conditions can be 
satisfied is illustrated by the results for a Timoshenko 
beam with linear and torsional springs attached at the ends. 
It is also illustrated that the inclusion of Timoshenko 
effects does not change the order of the interpolating 
function, which has been done in almost all the earlier 
formulations. The rotational stresses pose no additional 
complexity when included in rings or discs. 

Finally, the conclusions are given in Chapter VII. 



CHAPTER II 


THE FINITE ELEMENT METHOD 

The general kind of the eigen-value problem, that is 
to be solved, can be formulated as follows. Consider some 
domain Dm bounded by the surface £ . Let 9 be a scalar 
function defined in the interior of Dm such that the 
behaviour of 9 in Dm is given by 

- f = 0 (2.1) 

where f is a known scalar function of independent variables, 
and 0^ is a linear or nonlinear differential operator, with 
physical parameters as known constants or known functions. 

2.1 METHODS FOR SOLVING THE PROBLEM 

The general problem is to find the unknown function 9 
which satisfies equation (2,1) and the associated boundary 
conditions specified on the surface J , 

The various method available can be classified into 
the following broad categories : 

i) Exact solutions 

ii) Approximate solutions. 

2.1.1 Exact Solutions 

Some of the techniques used for finding the exact 


solutions are 
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a) Separation of variables 

b) Similarity solutions 

c) Fourier and Laplace transformations. 

The exact solutions are possible for a small 
number of relatively simple problems. Most of these have 
already been solved. They are the classical problems, 

2,1.2 Approximate Solutions 

The majority of the engineering problems are solved 
by the approximate solution procedures. Some of these 
procedures are s 

a) Perturbation methods 

b) Power series 

c) Probability schemes 

d) ’ Method of weighted residuals 

e) Finite difference techniques 

f ) Rayleigh-Ritz method 

g ) Finite element method. 

The regular and singular perturbation methods are 
applicable only when the nonlinear terms are small compared 
to linear terms. The power series method is powerful, but 
the generation of coefficient for each term is a tedious 
work. Also, at times it is difficult to show the convergence 
of tlie series. Hence, the usefulness of these methods is 
limited. The probability schemes are useful only when the 
quantity desired is a statistical parameter. 
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With the use of high-speed digital computers the four 
currently outstanding method for obtaining approximate solu- 
tions of high accuracy are the method of weighted residuals, 
the finite difference method, the Rayleigh-Ritz method, and 
the finite element method. The approximating functions, to 
be used in the method of weighted residuals, and the 
Rayleigh-Ritz method, must satisfy the boundary conditions 
on X , The pre-satisfaction of the boundary conditions 
makes the selection of the approximating function awkward for 
complicated/stringent boundary conditions. Hence, the 
use of these methods is again limited to the problems with 
simple boundary conditions. 

In the finite difference and the finite element 
methods, the selection of the function is not governed by 
the boundary conditions to be satisfied. Instead, the 
conditions are applied later to the problem. The appro- 
ximation of the complicated geometries by the finite 
difference method is more crude as compared to the finite 
element method, because in finite difference method the 
domain is sub-divided with the help of parallel straight 
lines only. Thus, the finite element method appears to be 
the most suitable approximate method out of the available 
lot. 

2.2 THE FINITE ELEMENT METHOD 

In the finite element method there are four ways in 
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which the elemental properties can be formulated. 

i) The direct approach, based on the direct stiffness 
method of structural analysis. The use of this 
approach is very limited because it can be applied 
only to relatively simple problems, 

ii) The variational approach is a more versatile and more 
advanced approach. It is based on the calculus of 
variations and involves extremizing a functional. The 
functional, in solid mechanics problems, is generally, 
the energy formulation in one form or the other, 

iii) A third and even more versatile approach to deriving 
the elemental properties is known as the weighted 
residuals approach. This approach begins with the 
governing differential equations of the problem and 
proceeds without relying on a functional of the 
variational statement. This approach is used for the 
solution of all the problems in the present work and 
is discussed in detail in Section 2,2.2, 

iv) The energy balance approach, based on the balance of 
thermal and/or mechanical energy of a system, requires 
no variational statements, 

2.2.1 General Procedure 

Irrespective of the approach used to find the elemental 
- properties, the finite element method always follows an 
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orderly step-by-step process. The various steps are : 

i) DISCRETIZE THE CONTINUUM : The first step is to 
divide the continuum or the solution region into the 
sub-regions, i.e,, the finite elements, 

ii) SELECT INTERPOLATION FUNCTIONS : The next step is to 
assign nodes to each element and then choose the inter- 
polation function to represent the field variable over 
the element. Quite often, the polynomials are selected 
as the interpolation functions because of the ease of 
differentiation and integration. The number of the 
nodes assigned to the element, the nature and number 

of the unknowns at each node, determined by the require- 
ments listed in Section 2.2.3, govern the degree of 
the polynomials. 

iii) FIND THE ELEMENT PROPERTIES : Once the finite element 
model has been established, the matrix equation expre- 
ssing the properties of the individual elements can be 
determined by the use of one of the four approaches 
discussed above, depending upon the nature of the 
problem. 

iv) ASSEMBLE THE ELEMENT PROPERTIES TO OBTAIN THE SYSTEM 
EQUATIONS : The next step is to combine the matrix 
equations expressing the behaviour of the elements and 
form the matrix equations expressing the behaviour of 
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the entire solution regionor the system. The nature 
and form of the system equations is same as that of 
individual element equations except that they contain 
many more terms because all the nodes are included. 

Before the system equations are ready for solution, they 
must be modified to account for the boundary conditions, 

v) SOLVE THE SYSTEM EQUATIONS t The assembly process of 
step (iv) gives a set of simultaneous equations, to be 
solved to obtain the unknown nodal values of the field 
variable, A number of standard techniques exist for 
solving linear as well as nonlinear simultaneous 
equations. 

Vi) MAKE ADDITIONAL COMPUTATIONS IF DESIRED ; Sometimes, 

the variables dependent on field variable are required. 
Now these calculations can be done with the values of 
field variables obtained at the nodes. 

2.2,2 The Weighted Residuals Approach 

Because of the broad choice of the weighting functions 
for error distribution principles, there is a variety of 
weighted residual techniques. The error distribution 
principle most often used to derive finite element equations 
is the Galerkin method. 

According to the Galerkin method the weighting functions 
are chosen to be the same as the approximating functions used 
to represent 9 . 
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Since equation (2.1) is valid over the entire 
solution domain, so it is valid over the domain of any 
element, e, also. Approximating the dependent variable, 9 , 
over the element as 


m 

9 “ 9^®^= I 9^®^ = LN^®^J{ 9 }^®^ (2.2) 

i=l ^ , 

where N^®^ are the assumed interpolation functions defined 
over the element and 9 ^®^ are the m-unknown nodal parameters. 
Substitution of equation (2,2) in equation (2,1) gives the 
residue as 


UP f ( e ) \ 

6-0(9 ) 


(2.3) 


Then, from the Galerkin method, we have, the equations 
governing the behaviour of an element as 

f d Dm'®' = 0 1=1,2,...,) 

f e "I 

(2.4) 


Since the operator contains higher order derivatives, the 
order of the interpolation polynomials required will be 
higher. To reduce the order of the polynomials we carry out 
the integration by parts in equation (2,4) till there is an 
overall reduction in the order of derivative. When integra- 
tion by parts is possible, it offers a convenient way to 
introduce the natural boundary conditions that must be satis- 
fied on some portion of the boundary. 
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2.2.3 Requirements for Interpolation Functions 

To guarantee the mono tonic convergence, the inter- 
polation functions in equation (2.2) be chosen so 

that the following general requirements are met, 

Huebner [54] : 

i) At element interface (boundaries) the field variable 
9 and any of its partial derivatives upto one order 
less than the highest-order derivative appearing 
under the integral sign in the modified form of 
equation (2,4), i.e,, after integration by parts, must 
be continuous. 

ii) All uniform states of 9 and its partial derivatives 
upto the highest order appearing in the modified form 
of equation (2.4) should have representation in 9^®^ 
when, in the limit, the element size shrinks to zero. 

The first one is known as the compatibility requirement, 
and the second as the completeness requirement , 



CHAPTER III 


BEAMS s FORMULATION AND SHAPE FUNCTIONS 

The beam is one of the basic components which exists 
in almost all structures. The effects of rotatory inertia 
and the shear were first studied in this component. 

Consider a beam of flexural rigidity El, modulus 
of rigidity G, cross-sectional area A, mass density P , 
and numerical shape factor for cross-section a. The coupled 
equations of motion for the deflection, w, and the bending 
slope, ^1/ , are [2,55] 




(3,1) 


P A 


3^w 

3t^ 


3X 


^ (aAG(|^ - t|<)) = 0 


(3.2) 


In case of a uniform beam these can be combined to give a 
single equation in deflection, w, as 


El 


8 2 , 




ft4 

w / .Pn *^W 


A pI+EI ;^) — + pi 

3t^ 3X^ 3t^ 


P 8^w 
“2 3t4 


= O 


(3.3) 


Using those equations, the finite element equations are 
derived and then the shape functions are deduced. 
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3.1 FORMULATION 

Uniform beam is considered first in Approach A and 
then the general beam in Approach B. 

3.1.1 Approach A s One Governing Differential Equation 
Assuming a polynomial function for the deflection as 


=l NKw) t"®' 


(3.4) 


over an element 'e'. Figure 3.1, and substituting in 
equation (3.3), we get the residue as 


Res 


= El -JV- + PA— (PI+EI -t) ^ g - 

3t^ 3x 3t^ 


3x 


+ pi 


p 


3t^ 


(3.5) 


Applying the Galerkin method to minimize the residue 
over the domain of the element, e, we get 


h 

/ 


N^ Res^®^ dx 


0 


i = 1 , 2 , , m 


(3.6) 


or 


El 


h 

/ 


N. T“ dx+ pA f N. 55 — 

J ^ 3t 


dx - 


3 X 








dx = 0 


X — f ^ f 


(3.7) 
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Integration by parts results in 


a3(e).h >2(e) 

El -(N| El ^ 


3x' 


9x" 


) 


h h 

+ El 


f NP — dx + 

i ^ 3 X 


e ) 


pA 


/ N. dx - (N.(PI+EI 


3x3t 


o 


( PI+EI / N! ^ 

O 


dx+ Pl ^ dx = O*, 

o 


i = 1,2, , , , , m 


(3.8) 


where ( ’) prime denotes a derivative w.r.to x and (*) 
dot denotes a derivative w.r. to time, t. 

Expressing equation (3.8) in vector notation and using 
equation (3.4), the elemental equations of motion are obtained 
as 


3 ( e ) 

k tA/'* ' 


({N} El AW-) 


3 X' 


El -t) 
aG'' 


2 (e) 




-({N‘} El ^ 

3x 


- ({N}(pl + 


o 


3 ( e ) , h ! \ 

"' ^) + El / (N") LN"J dx + 

I o 

o 


3x 3t' 


h 


‘aG 


o 


or 


h 

+ EI^) /(N'} 

LN’ Jdx)|w/^®^ 

0 


(ne) ^ Q 

(3.9) 



(3.10) 
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where 

= EI[A^] 

[M^]<®>= pA[A3] + ( PI+EI [Aj] 

[M2 ](®>= pi ^ [A 3] 
h 

[Aj^] = / {N*' } LN"J dx 

o 
h 


IA 2 I = / [{N«} [N*J dx 


o 

h 


[A3] = / {N } LNJ dx 


and 




( ne). 


9 V") 


-({N} El 


3je) 


■w'”' ih 


+({N }(pI+EI 1-^) 


3x 3t" 


+ 


o 


+({N»} El ) = {Y. } + {Y. > 

3 x^ 'o ^ 


3.1#2 Approach B : Two Simultaneous Governing Differential 
Equations 

In case of a non-uniform beam, the coupled equations of 
motion (3»l) and (3.2) could not be combined into one equa- 
tion. An alternative approach to solve these equations 
simultaneously is used. 
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Assuming a polynomial interpolating function for the 
deflection, w, as 

= Ln”j (3.11) 

and another polynomial interpolating function for the 
bending slope, , as 

= LNh (3.12) 

where the superscripts w and denote the interpolating 
functions for and respectively, over the 

domain of the element *e’. Substituting equations (3.11) 
and (3,12) in equations (3,1) and (3.2), the two residues 
are obtained as 


± (El 


pi- 


.2ie) 




(3.13) 


(3.14) 


The residue Res^®^j equation (3.13), is from the 

equation of motion of rotation, therefore, it has to be 

minimized with respect to the interpolating functions of the 

( e ) . 

bending slope. The residue, Res^ of equation of trans- 
latory motion, equation (3.14), has to be minimized with 
respect to the interpolating functions of the deflection. 
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Hence, applying the Galerkin method, we get 


/ dx = 0 i = 1,2, 1 (3.15) 


o 

h 


/ nV Res^®^ dx = 0 i = 1,2, ..., m (3.16) 


o 


where 1 and m are the number of nodal parameters for bending 
slope and deflection respectively. 

Substituting equation (3.13) in equation (3,15) and inte- 
grating by parts, we get 

^ h (e) ^ 

f PT M dv 4 - f K. 

1 


El 


- / NI El 1^4^ ^ 

o o ^ ^ o ^ 


3X 


~ / N. PI Tj; ^ ^ dx = 0 j i = 1,2, .. .,m 

O 

(3.17) 

Expressing in vector notation and using equations (3.11) and 
(3,12), the above equation (3.17) reduces to 

h 


({N> - / EI{N’} LN’ J dx {4»} - 


3 x 


h 


0 o 


n 

/ aAG{N} in J dx{t{;}^’^®^ + / aAG{N}LNj dx {w}^’^®^ 

o 

^ = 0 


/ pi {N > LNJ dx y 


(3.18) 


0 
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Similarly, equation (3.16) yields 

a AG + / aAGlN’'^ }LN’^Jdx - 

^ |o o 

h h 

/ aAG {N''"}lN J dx + / PA Cn'^ } [N'^Jdx = 

o o 

(3.19) 


Equations (3.18) and (3,19) can be combined as 


— 

[M^] 

[ 0 ] 

(e) 

1 l(ne) 
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! 

r- 
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(ne) 
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(3.20) 


where 


h 


[K,] = 


/ aAG {N' } LN’ J dx 

o 


h 


[Kp] = - / aAG {N’'" } dx 


[K3] 


= / (El {N* LN> j + aAG { n’^} [n’^) dx 


o 

h 

[M^] = / pA {n'^ >[n'^J dx 

o 

h 

[M^] = / P I { N } LN dx 


o 
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r 


and 


(ne) 
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({N^ } aAG( 
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(e) 


9x 


1 


o 


,(e) ih 




J 


3.2 SHAPE FUNCTIONS 


A review of the various shape functions reported 
in literature, as classified in Table 1.1, gives a better 
appraisal of the shape functions proposed in the present 
work. 

i) Element A [5-ll] s As reported by Thomas et al [4] 
that the element A formulations [5-ll] are equivalent and 
can be reduced to anyone of them, so as a representative, the 
derivation given by Przemieniecki [ll] is taken. The mass 
matrix has been derived using the kinetic energy formulation. 
The equations of static equilibrium have been used to 
derive the stiffness matrix by what is known as the direct 
method. This results in an inconsistent stiffness matrix, 
however, the error introduced is small enough to effect the 
accuracy of the results significantly. 

ii) Element B [l2,13] i Though the nodal parameters used, 
in this case, are w and ^ but the elemental matrices trans- 
form to those given by element A, with, the help of suitable 
coordinate transformation. 
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iii) Element C [l4] s This element has two nodes with 
eight nodal parameters. It has a limitation that the 
abrupt changes in cross-section cannot be taken into 
account. There is a discontinuity of shear at such 
sections, but the element enforces continuity. Moreover, 
a higher order of interpolating function for shear has, 
unnecessarily, been employed. 

iv) Element D [l3] : In this case three nodes with 
seven nodal parameters are used, A quadratic interpolation 
for shear is used. This element has the same drawbacks 

as those of element C. 

v) Element E [l5] : The cubic interpolation functions 
for both w and i|» have been used. There are four nodes 
and eight parameters. Since the continuity of shear is 

not enforced, this element avoids the difficulty encountered 
by elements C and D. 

vi) Element F [4] t It has two nodes with w, and the 
shear deformation, y> as the parameters. Using the relation 
[2,55] 

If = ♦ + Y (3.21) 

with the help of suitable transformations, it can be shown 
that the parameters reduce to w, ^ and y> i.e.» a cubic 
interpolation for the deflection and a linear for shear. 

This element also has the same difficulty as the elements 
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C and D, namely, this also enforces the continuity of shear 
at the abrupt changes of the cross-section. 

A careful study of the above elements reveals that 
excepting elements A and B, which have inconsistent matrices, 
all others have been derived using interpolating functions 
of higher order than that required by the problem. More- 
over, no reason has been given for selecting the particular 
functions. Each element has one or the other limitation, 
hence, none of them is general enough to solve the problem 
in its entity. 

In the present analysis the shape functions are selected 
by considering the compatibility and completeness require- 
ments, Huebner [54], 

3,2.1 Shape Function for Approach A 


The highest order of derivative under the integral 
sign in equation (3.9), is two, hence, the compatibility of 
the variable w and its first derivative is required. The 
completeness of all derivatives upto the third order, the 
highest order appearing in equation (3.9), is required. 
These are same as the requirements for standard interpola- 


ting function for the Bernoulli-Euler beam. Thus the fun 


ction is, Desai and Abel [56], 
= [C^(3-2^) 
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where 

^ = 1-x/h and ^2 = (3.23) 

The matrices [A^] , [A2] » [A^] and the vectors (Y^ ), {Y2 ) » 
required in equation (3.10) using this function, equation 
(3.22), are given in Appendix A. 

3.2.2 Shape Functions for Approach B 

The highest order of derivative of w and under 
the integral sign in equations (3.18) and (3.19) is one, 
so compatibility of w and iff is required. The completeness 
of w, and their first order derivatives, the highest 

order in equation (3.18) and (3.19), is required. The 
functions satisfying these are 


= L ( 2 ^ 

VkJ. 

( 3.24) 

= L 5^ SjJ 

i 

> 

1 

( 3.25) 

Therefore, 



LN’^'j = LN^^J = LNj = 

L 5^ C2J 

(3.26) 


The matrices [K^] , [K^]? » [%] » and the 

vector {82/’^®^ using these interpolating functions and 
for a linear variation of cross-sectional properties, are 


given in Appendix A 



33 

With the present formulation any variation of cross- 
sectional properties can be handled, without any approxima- 
tion, at the time of integration of the elemental matrices. 
This is a clear-cut advantage over the earlier formulations 
where assumptions have been made to incorporate the cross- 
sectional property changes. 



CHAPTER IV 


RINGS AND SECTORS j FORfvlULATION AND SHAPE 

FUNCTIONS 

4.1 EQUATIONS OF MOTION 

The equation of motion of a thin ring, without 
rotational speed, has been given by various authors [22,57], 
To include the effects of rotational speed, n , in the 
equation of motion, consider an element m m^ n^ n of a 
rotating ring or sector, as shown in Figure 4.1, of 
radius R, flexural rigidity about first axis EI^, cross- 
sectional area A, and mass density p , If T is the 
tension, F is the shear force, M is the bending moment, 
q is the uniform load acting on the element, u is the 
radial displacement and w is the tangential displacement 
then the equations of motion in tangential and normal 
directions are 

- F = pAR (4.1) 

3t 

( F+ d-0)cos — ^ — F cos — ^(T+ d©) sin — ^ + T sin — ^ - 

q( R d©) = p AR[ " ‘ ^^ ' + (R+u) flP ] d© 

3t^ 


(4.2) 
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and, the moment about mass centre gives 


36 


+ FR = 0 


(4.3) 


Simplification and taking limit as d6 o, reduces equation 
(4.2) to 


ap 

30 


+ T - qR 


PAR[— & 
3t^ 


+ (R+u) 


(4.4) 


From theory of elasticity, Prescott [58], we have 



(4.5) 


and condition of inextensibility of centre-line gives [22] 


u 



(4.6) 


Using equations (4.5) and (4.6), equations (4.1), 
(4.4) can be combined as 

6 4 2 4 2 2 

a°w . ^ a w . 3 w PAR^ r t a w 
^ + 2 J + — TT — - w; - 

ae° 30^ ae ei^ 3t 

ia oAR^ o2 
Ell 36^ 


(4.3) and 


(4.7) 


For free in-plane vibrations of an inextensible, rotating 
thin ring or sector equation (3,7) reduces to 
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39 ° 3 0 ^ ^-*-1 36 ^ 


4 2 

PAR^ 3 


El 


2 

/ 3 w 


1 3t^ 30 


w) 


(4.8) 

which is the governing differential equation of motion. 


In the governing differential equation of motion for 
a non-rotating ring or sector including the Timoshenko 
effects [22], the rotational speed effects can be included 
in a similar fashion as 
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3t^ AR^ 


) ii" + 

304 


(1- ii-Ca- Ig)) - (1+ l37)w)) (4.9) 

AR^ 36'^ AR"^ 

4.2 FORMULATION 

A typical finite element of a ring or sector subtending 
an angle p at the centre, is shown in Figure 4.2* The 
coordinate system is as shown and, hereafter, © will refer 
to local coordinate unless otherwise stated. 

Assume a polynomial function for the tangential 
displacement as 

= [NJ (4.10) 
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4.2.1 A Thin Ring or Sector 

Substituting equation (4,10) in equation (4,8), we 
get the residue as 
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where 


and 


X = 1+ 
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Applying the Galerkin-method to minimize the residue 
of equation (4.11) over the domain of the element, we get 
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Integration by parts results in 
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(4.13) 


d© = 0, i = 1,2, ..., m 

where ( ' ) prime denotes a derivative w.r. to 0 

and ( *) dot denotes a derivative w.r, to time, t. 

Substituting equation (4,10) and using vector notation in 
equation (4,13) we get 
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+ [K]^®^ (4,15) 
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where 

= C([A3] + [A^]) 

[K](®^ = [Aj^] - 2[A^] + XEAj] 

such that o 

P 

[A^] = / {N»” }LN’” j de 

o 

r n 

[A2] = j {NMLW'J de 

o 

p 

[A3] = / {N’ }LN’J de 

o 

p 

[A4] = / {N } LNJ de 

o 



= {Y^> - {Y2 > + {Y3 } 

4.2.2 A Thick Ring or Sector 

Substituting equation (4.10) in equation (4.9) we get 


the residue as 
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3 


.4je) ^2je) 

- + X 


3 e 




3^" , 3-w 


r(- 


,2 (e) 


4^ o 
3t^ 30^ 


’■ 30" 

(e)^ 3^ 

w^ 0+ — 5 
3t^ 


AR' 


1 ( 1+ i_) 


(1- ^(2- |g)) -£4^ - (1+ 

AR 30 AR 


3© 

)) (4.16) 


Applying the Galerkin method to minimize the residue of 
equation (4.16) over the domain of the element ’e’, we get 




/ N. Res^®^ de = 0 i = 1,2, 


• « • 9 


m 


(4.17) 


or 


/ Ni(-^ 


6,.,(e) A,(e) ,2,.,(e) ,e ,h ,2 ^ s'* , sV ® > 


W' 


o a " , V a ^ J / 1 r J / 3 ^ 

30^ 3 0"^ 3 0^ A? 3 0^ 


~~( _ — : 

at^ AR‘ 




36' 


(e)\ . ^ ( '*'1 i' H ■ E N "^^l E ^^ 3 w 


,2je) 


^2- Ig^) ^ 

AR^ 3 © 


11 4®) 

Integration by parts results in 


( 14- "7? ) w )) d© = 0y i — 1^2, ...yin 

AR' 



43 


(N^(' 






I p/E ( 1 

AR 


2 


(Nl(^ 

^ 36^ 


4. 


2 ^ ^ 
302 


II 

C -%! 
AR 


3 

aG g0^ 


^ se'^ 


- ;(Nr^ 

3 e'^ lo d ^ 3 


^ - 2K': 


■^^4“ + X n! M^)de 

» q 2 13© 


3 e' 


c^ / (N^ 


S-(X)' 

aG'^T^^ 

AR 


+ w 


w ^ ® ^ ) d© - 



Substituting equation (4.10) and using vector notation in 
equation (4. 18), we get 
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> 3 je) 


+ 2 + X 


x^ 3 'w(e) 

^0 + 1 ©~ 


3;A(e) 


X-f 1 + 3 w^ 

AR 5 0'^ 




({N’ .}( -5.^!^- + 2 C -4 

a@^ a e^ ar 


I>, C n 2 **( G ) I B 

( 1 + |g) — V ‘ ) ) + 

AR 3 ©^ o 


-3 (e) p p 

({N«* } -^-\ - ) - f ({N*»}LN"‘ J-2 {N" }LN»J + 

o o 

I P 

X{N’} L»' J)d0 {w } - |pr(-^)^ / (<N' } LN' J + 

AR o 


p I 

{ N > LNJ)d0 ” ^ / (“^X+ Iq) {N» } LN"J 

/ £ I / £ \\. o . , .AH. . 


svo ,y-in 

))l 


J ■+■ 


( 1+ -^) {N } LNJ)d@ = 0 

AR^ 


(4.19) 


+ [ k ]^®3 { w /'"®3 = {82}^’^^^ 

(4.20) 


where 


[M,] 


2 ^'®^ = is c 2 ([A 3 ] + [A^]) 


= C(-^ 1+ |g) [Ag] + (1- -3(2- iG))[A3]+(l+ ^)[A4]) 


[K]<®) = [A^] - 2 [A 2 ] H- XCAj] 
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and 


{^2} 


(ne) _ 


(W }(i!42^ + 2 1^4^ + X + 

3 6^ 8 ©'^ 


3© 


aG 


3W-(-) ^ 

70“ 




Ar2^^'^ aG^ ggS 


3,7,(e) 


A? ^ 1 

^ O 


- ({N-} + 

3© 


3^w^ ® ^ 
1?~ 


C ( 1+ 
AR 


E ^ 3^#®^ 

aG^ ■;i5— 


)) 


+ ({N"}- 


3 V® \ jp 


3©' 


= {y^ } - { Y^ > + { Yr } 

4.3 SHAPE FUNCTIONS 

The shape functions for analysis of thin rings and 
sector reported by various authors have been compared by 
Sabir and Ashv^ell [34]. The shape functions are 

i) A two-dimensional (ring) form of Cantin and Clough’s 
shape function for cylindrical shells [27] is 


u 


w 




2 3 

sin© + a 2 cos© + a^ R cos ^-p sin© + a^y + a^y 

(4.21) 

cos© - 32 sin© - a 2 R( l-cos cos©) + a^y 

(4.22) 

(4.23) 


where y = R© 
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ii) A two-dimensional form of 48-nodal parameter shape 

function of Bogner et al [29] for cylindrical shells 


iii) 


u = ai + ajY + ajY^ + 

2 3 

w = a^ + a^Y + a^y + a^y 

has 8-nodal parameters, namely, w,w' - u/r, 
at each node. 

A reduced form of the above shape function, 
(4.24) and (4.25) is 

u = a^ + 327 + a^y'^ + a^y^ 
w = + 357 


(4.24) 

(4.25) 
u and u' 


equations 


(4.26) 

(4.27) 


having 6 nodal parameters, 

and 

iv) A constant circumferential strain, e, and linear 
change of curvature is given by the model [33] 


u = a^ cos© + 32 sin© + a^ - a^^© (4.28) 

2 

w = -a^sin© + a2 cos© + 82 + a^Q + ^a^Q (4.29) 

It is obvious from equations (4,21) through (4.29) 
that none of the shape function (i) to (iv) satisfies the 
equation of inextensibility of the centre-line of the ring, 
equation (4.6). In fact, the number (iv) shape function 
is based on the assumption of linear circumferential strain. 
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Thus none of these shape functions represents the problem 
under consideration, which results in convergence problems, 
[34]. 

In the present work the shape functions are deduced 
considering the compatibility and completeness requirements 
[ 54 ] as governed by the equations (4.14) and (4,19). The 
highest order of the derivative appearing under the inte- 
gral sign in both the equations (4.14) and (4.19) is three, 

( o ) 

thus, a compatibility of all derivatives of w upto 

second order is required at the element interface. The 

) 

completeness of all derivatives upto the fifth order, the 
highest order occurring in equations (4.14) and (4,19), is 
required. Since the requirements for the interpolation 
functions for the analysis of thin and thick rings is same, 
the functions must also be the same. A function satisfying 
these requirements is 

= a + b© + c@^ + d6^ + + f©^ (4,30) 


Expressing in non-dimensional parameters, and 3^*^ 

nodal parameters, the tangential displacement function becomes 


J- ® c3( 10-15 Cj^+6 

10 - 15 ^+ 6 ^) 


-P ClfcC i+3ei) 


w. 

J 

I 


w. 


w 


u 


k 

I 

k 

It 


k 

(4* 31) 
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where = 1-e/p and ?2 = ®/p» 

The various,;, matrices and natural boundary condition 
vectors, evaluated using tho'‘'above shape function, equation 
(4.31), are given in Appendix B. 



CHAPTER V 


DISCS : FORMULATION AND SHAPE FUNCTIONS 

The equation of motion for free vibrations of a 
rotating disc of thickness T, including Timoshenko effects 
is [58-60] 


D v\-{pl4-D pi 




+ pT = 


1 


(5.1) 


where 


D = El/(l- v^)» the flexural rigidity 

W=W(R,e,t) 

cr©Q are the stresses induced due to rotation. 

All other symbols have the same meaning as in Chapter III. 

Assuming a cosine function for the variation of displace- 
ment in the angular direction, we can write 


W = w(R,t) cos n0 


(5.2) 


Substituting equation (5.2) in equation (5.1), we get 


D vV(pI+D v2k+ pi ^ w- + PT w = I ^(RP + 


n TOgg w 


(5.3) 
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as the one-dimensional equation to be solved for free 
vibrations of a disc. 

5.1 FORMULATION 

A typical annular element, e , of a disc is shown in Fig- 
ure 5.1. The coordinate system is as shown. Transforming 
equation (5.3) in terms of the local variable, r, we get 
after suitable algebraic manipulations, 

^(^2((R «) -^^)-Cl+2n2) 2^42! ). 

j 3^ (R.+r)'^ 


•J J 

^ ” R^r = 0 

R^+r ar j rr ar (d ^)2 


(5-4) 


Assuming a polynomial function for the displacement 


= [NJ 


(5.5) 


and substituting in equation (5.4) we get the residue as 


Rej 




.2 (e) 


(e) 


■j ~ 9r ar U 

-4^-2 D -i -vC®) n^w^®^ 


w(®h-(PH-D :4)-5^j|;:((R,+r) |S 


( Rj+^)' 


aG^Rj+r^B r^ ^ j ^ ar 


R .-t-r 

sj 


)4- 


pl aG P 


i^^TcTqqW 
(Rj+r)' 


P (e) 
(e) 


w + oTw^®^'- ~r ((R,-+r)cr ^ )- 

Rj 4 -r ar ' ' 3 rr ar ' 


(5.6) 
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Applying the Galerkin method to minimize the residue 
of equation (5,6) over the domain of the element, we get 


h 2% 

/ / N. Res^®^(R.+r)d€) dr = 0 


= 1,2, m (5.7) 


o o 


/ N.(R.+r) Res^®^ dr = 0 i = 1,2, 

X J 


(5.8) 


Substituting equation (5.6) in equation (5.8), we get 


f I? ® ^ 2 3 1 ® ^ 

D / N^(7-^(F^i+r) — 5-)-(l+2n ) ) + 

<3 


1' . Z' " 1 

o 3r 


O W^®^)dr - ( Pl+D^) / N,(|r((R,+r) 


(R.+r)- 


^ aG^ ' 3r 

2 "(e) oh t\ 

)dr + PI ~ / (R,+r)N/W'^®^dr+ 


Rj+r 


aG ' 


pT / (Rj+r)N^w^®^dr - T / ^ ^1+^^^rr If " 


i 3r" " j '"rr 3r 


n^T / N.w^®^ dr = 0 i = 1,2, .,., m 

o 


(S3) 


Integrating by parts and using equation (5.5) and vector 
notation, we get 



such that 


h 


[A^] = / (R.+r){N” >LN” J dr 


o 

h 




[AJ 


o 3 
h 

f Tr^sW} LNJ dr 


[A4] = / ar^^(Rj+r) {N’ }LN'J dr 


o 

h 


[Aj = / 


^ee 

Rj+r 


{N}LNJ dr 


h 

[A^] = I (Rj+r) {N» } LN*J dr 
o 

h 

[A7] = / (R-+r) {N} LNJ dr 

[Ag] = / R^{N)tNJ dr 

r\ 


and 

B 


(dd) ^ .(jN j(D(±{(Rj+r) 

- (pI+D ff" )) 


R 


T(Rj+r)cr^^ 3r 


32je) 


({N» }D(R.+r) 

J ar 


= - {Y^ } + {^2 } 
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5.2 SHAPE FUNCTIONS 

The simplest of the shape functions used by the 
various authors [43—49] to study the effects of rotation 
and pre-stress on free vibrations of thin circular discs, 
is, the transverse displacement, w, is approximated by 

w(r,e) = (a^+a 2 r+a 3 r^+a 4 r^)cos nO (5.12) 

and the radial displacement function, u, is 

u(r) = b^+b 2 r+b 3 r^+b 4 r^ (5.13) 

Since cubic expressions are used for the radial and trans- 
verse displacement, a total of eight nodal parameters are 
required per element. There are other shape functions 
which require still higher number of nodal parameters per 
element. 

The shape functions for the analysis of a thick disc, 
proposed by Wilson and Kirkhope [50] is 

2 2 

w(r,©) = (a^+a 2 r+a 3 r 4-a4r )cos n© (5.14) 

Yp(3:,©)= (a^+a^r) cos n© (^.15) 

and 

(ay+agr) sin n6 (5.16) 

In the earlier works, reasons for selecting the nodal 
parameters or the shape functions are not given and in almost 
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every case the number of nodal parameters or the order of 
the functions is selected, unnecessarily, more than the 
required. In other words, the functions have been 
selected arbitrarily. 

The compatibility and completeness requirements for 
the interpolating function [54] are governed by the 
equation (5.10). The highest order of derivative under the 
integral sign is two, so the compatibility of transverse 
displacement, w, and its first derivative is required. 

While the completeness of all derivatives upto third order, 
the highest order in equation (5.10), is required. These 
requirements are same as for the interpolation function for 
the beams, Chapter III. Hence, the function given by 
equation (3.22) is the required function, namely, 

Iw.’ i 

^ > 

w'®' = Ui( 3-253^) £352^ 4 ^ 3 - 2 ^ 2 '> r*' [ 

■v. ^ , 

The elemental matrices [A^] to [.ig] and the vectors {Y^ } 
and {Y 2 } using this function are given in Appendix C. 

The requirements for the interpolation function remain 
the same even if the Timoshenko effects are neglected. Hence, 
the function for both cases remains same. 



CHAPTER VI 


RESULTS iW) DISCUSSION 

Before discussing the numerical results, a brief 
discussion of the solution procedure is given. 

6.1 SOLUTION PROCEDURE 

An examination of the elemental equations (3.10), 
(3,20), (4,15), (4,20) and (5.11) reveals that the 
assembled system equatiors vvill essentially reduce to one 
of the following equations, after application of the 
boundary conditions. 

i 

6.1.1 Classical Equations 

V^hen the Timoshenko effects are not considered, the 
governing equations obtained are the simple second order 
equations of the type 

[M] + [K] ={0} (6.1) 

The dynamical matrix for such equations is well esta- 
blished and is constructed as [55,62] 

[D^] = [K]-^ [M] (6.2) 

6.1.2 Equations with Timoshenko Effects 

A fourth order governing equation of the type 

[Mg] fw}'") + [M^] + W = {0} (6.3) 
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is obtained, when the Timoshenko effects are included in 
the analysis. The dynamical matrix for the fourth order 
equation (6,3) is constructed by the direct method, 
Przemeineicki [ll], as 

[K]”^ [M^l [K]“^ [M2] 

[^ 2 ^ = ( 6 . 4 ) 

~[I] ' [0] 

In case of the structural components having rigid-body 
modes, a shift of origin, or the o 

Shift, Bathe and Wilson [61], is applied to the stiffness 
matrix. The shift is applied to avoid the singularities 
while inverting the [K] matrix, 

6.1,3 Solution Algorithm 

The method of matrix iteration is applied to the 
dynamical matrix for finding the eigen values and^eiven- ) 
vectors of the system of equations [55,62,63], 

In this method, we assume a starting vector and multiply 
it with the dynamical matrix. Compare the new vector with 
the old vector for convergence. Drop the old vector and 
repeat the process with new vector till the two vectors 
become almost equal. 

The solution procedure can be described by the following 


s teps 
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i) Get the assembled matrices from the elemental matrices 
and apply the appropriate boundary conditions, 

ii) Apply the Shift, if required, 

iii) Construct the dynamical matrix, [D] , from the 
assembled matrices. 

iv) Select an initial vector i and set i = 0, 

v) Construct a new vector as }, 

vi) Compare with {X^}, <e for all 

k, stop, otherwise, set i i+1 and repeat from 

step (v) onv^ards. 

To save the computational efforts, the symmetry was 
exploited wherev.c'r possible. The free-ring results were 
obtained by analysis of a 180° - sector. The constraints 
on radial displacement were applied from the symmetry consi- 
derations. This also avoided the problem of singular stiff- 
ness-matrix. 


6.2 NUMERICAL RESULTS AND DISCUSSION 


The various numerical results obtained are discussed 
in this section. The follovdng non-dimensional parameters 
are used throughout the discussion to make the results more 
general* 


•sr^ _£al:^2 
n " EI n 


2 I 
re = 

klT 


(6.5) 


( 6 . 6 ) 
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se 


2 


£I 

aAGL"" 


(6.7) 



( 6 . 8 ) 


k 


t 



(6.9) 


_ PAL^ ^2 
“ £I 


( 6 . 10 ) 


The results of various authors have been reduced to these 
non-dimensional parameters for comparison. The results 
for many of the unattempted problems are also obtained, 

6,2.1 Beams j Approach A 


A number of problems have been solved and the compari- 
son of results has been made with earlier results wherever 
possible. 

Table 6,1 shows the results for the cantilever beam 1 
analysed by Thomas ct al [4], The non-dimensional parameters 
for the beam are re = 0,002 and se = 0,007, It can be 
observed that the convergence rate is quite fast. 

In Table 6,2 a comparison of error in results has 
been made with Thomas et al [4], for the beam 1 [4], Mag- 
nitude of 9^ error in the results of Thomas et al [4] has 
been read from the graphs. 
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Table 6.1 

Natural frequencies,"^^, for cantilever beam 1 
of Thomas et al [4] 


Mode , n 

Exact results 
[3] 

Finite element 

results 


2 elements 

4 elements 

6 elements 

1 

3.5155 

3.5173 

3.5157 

3.5156 

2 

22.0155 

22.2020 

22.0412 

22.0210 

3 

61.5712 

74.9747 

62.0465 

61.6837 

4 

120.4468 

216.2965 

122.1886 

121.2162 


Table 6.2 

Comparison of error in results 

and present work 

of Thomas et al[4] 

Mode,n 

error in results of 


reference[ 4] 

present 

re fere nee [4] 

present work 

9 nodal para- 

work, 8 

18 nodal 

18 nodal ' 

meters 

nodal para- 
meters 

parameters 

parameters 

1 0.01 

0.006 

~ 0.00 

0.000 

2 0.30 

0.100 

0.05 

0.006 

3 1.25 

0.800 

0.20 

0.040 
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A better rate of convergence is observed from the 
Table 6.2 for the present formulation as compared to the 
earlier formulation, element F, claimed as the best [4]. 

The results for the cantilever beam. 2 of Thomas et al 
[4] are presented in Table 6*3, the non-dimensional param.eters 
are re = 0.049 and se = 0.100. A good rate of convergence 
is observed even for the higher values of re and se. 

Table 6,3 

Natural frequencies ,”15^, of cantilever beam 2 of 

Thomas et al [4] 


Mode,n 

Exact 

results 

[3] 

Finite 

element results 


2 elements 

4 elements 

6 elements 

1 

3.4194 

3.4204 

3.4191 

3.4190 

2 

18.6336 

18.7693 

18.6668 

18.6545 

3 

44.6955 

52.6259 

45.1671 

44.9597 

4 

74.6361 

107.2538 

76.8465 

75.9180 


A comparison of ^ error in results of the cantilever 
beam 2 of Thomas et al [4] with present work is made in 
Table 6.4, 

The present formulation shows a better rate of conver- 
gence and a higher accuracy of results with fewer nodal 
parameters compared to the earlier formulations. 
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Table 6.4 


■ Comparison of ^ error in results of reference [4] 

and present work 



,.,Mpde,n.- 

error 

■ 

in results of ?■ 



reference ' t 4] ' 

9 nodal parameters 

present work . ' 

. ; ; 8 nodal ' parameters ' 


i 1 • ■ .■ 

- 

■ 0^01 


2 

0.5 

0.’l9' ■* 


3 

.x6..25 

■ ■ k . ■ /’ ' ,, 

1.06 1 



The io error in the results of Kapuir [l4] for' cantilever 
■ and simply supported beams are compared with.- ^ error in 
> results of present work in Ta^le., 6.5. . The non-dimensional 
parameters art re = 0.080 and se =■. 0.160. 

There seems'.to be some error in Kapur's results [l4] 
for 4th mode of the simply-supported be.am, because^for the; 
higher rr?odos- .thd error increases mongtonically in ’all ; 
cases. The accuracy of. the present formulation is ■'in-general 
higher than the earlier ones, though the error is%lightly 
more in first two modes of the cantilever beam, 

■A' comparison of the ^ error in' results of present 
v/vork and Davis et al [lO] shows -Ui at there is a significant 
improvement in the accuracy. . Ee-sultte are presented in 
fdblG 6.6. In this case also the ^ error for the fourth mode 
of the simply-supported beam, as quoted by Davis et al [ 10] 
seems to be in error. 
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Table 6,5 

Comparison of io error, in Kapur’s results [14] and 
the present work, for 2 elements 


Mode,n 

i<> error in Cantilever 
beam results 

io error in Simply- 

supported beam results 


refercnco[l4] 

present v>/ork 
work 

re ference 

[14] present 
work 

1 

0,03 

0.0790 

0,29 

0.0266 

2 

0.57 

0.8387 

6.00 

0.4060 

3 

7.88 

, 3.0130' 

9.50 

1.8950 

4 

14.90 

7.3375 

8,00 

12.4750 


The natural frequencies of simply-supported (S-S) and 
clamped-clamped (C-C) beams for low values of the non- 
dimensional parameters, re =,0.02 and se = 0,02, are obtained, 
and are given in Table 6,7 along with the exact results, using 
frequency equations of Huang [3], A similar table. Table 6,8, 
for higher values of non-dimensional parameters, re = O.iO and 
se = 0.10, is also obtained. An excellent trend of convergence, 
for both higher and lower values of re and se, is observed, 
except for an isolated case of a fourth mode of clamped- 
clamped beam for higher re and se. This may be due to the 
ill-conditioning of the equations. 



Table 6.6 

Comparison of i> error in results of Davis et al [ 10] and the prosont work 
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Table 6.7 

Natural frecjuencies simply— supported and clamped-clamped beams 

re = 0.02 and se = 0.02 
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Table 6.8 

Natural frequencies,”**^^} of simply— supported and clamped-clamped beams, 

re = 0.10 and se = 0.10 
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The vers all ty /of the formulation is proved by 
obtaining the natural frequencies of a few cases of non- 
simple boundary conditions, e,g,,when different types of 
springs are attached to the free end of a cantilever beam. 
Figure 6.1 gives a plot of “/“* versus k. for the first 

o 

three modes of ,a cantilever beam when a linear spring, k^, 
is attached, where w* is the frequency for k =0, the 
value of non-dimensional parameters are re = 0,02 and 
se = 0,02, The additional spring has maximum influence on 
the fundamental mode and almost negligible on higher modes. 
An •:.;identical trend can be observed for re = 0.07 and 
se = 0.10 values, as shown in Figure 6,2, 

Figure 6.3 shows the variation of frequency when the 
same spring is attached to beams of different dimensions, 
such that se = 2 re, is the frequency for re = 0.0. 

The trend is same as for k = 0, except that for lov^er modes 

o 

the rate of decrease of frequency is not so steep. 

Figure 6,4 shows a similar graph when an additional 
torsional spring, k^, is attached. In this case the rate 
of decrease of frequencies has been reduced for all modes 
compared to that for k^ = 0. 

6.2,2 Beams s Approach B 

The natural frequencies of uniform cantilever (C-F) 
and simply— supported beams are presented in Tables 6.9 and 
6.10 for low and high values of non-dimensional parameters. 
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re and se. For low values of re and se the convergence rate 
is extremely poor. Table 6,9, whereas, the convergence rate 
is quite satisfactory for the higher values of re and se. 
Table 6,10. This poor rate of convergence is due to ill- 
conditioning or locking, as sometimes called in literature. 
Similar convergence problems were encountered in static 
analysis by various authors, Gallagher [64] and Hinton and 
Owen [65] , To avoid this ill-conditioning, the domain was 
divided into elements of unequal lengths to change the 
nature of the simultaneous equations to be solved. But 
there was no improvement. In static analysis the ill- 
conditioning was removed [65] by the reduced integration 
using l-point Gauss-Legendre quadrature formula [6l] for 
the shear stiffness matrix. The results using this 
technique for [K^] in equation (3.20) are also given in 
Tables 6.9 and 6,10, and the marked improvement in con- 
vergence rate can be observed. Still the effects of the 
ill-conditioning can bo observed for the fourth mode of 
the simply-supported beam, re = 0,10 and se = 0.15, where 

the frequency has converged to a value less than the exact 

t 

value. 

The natural frequencies, 1^, of tapered cantilever 
beams are given in Table 6.11, The taper ratio, Tp, as 
shown in Figure 6,5, is 0,90 and 0.70, re = 0,10, and 
se = 0.15. As expected the frequency increases compared to 
the uniform beam. A good rate of convergence is observed. 
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6.2.3 Thin Rings 

Non-dimensional natural frequencies, are obtained 

for various configurations of rings and sectors. Natural 
frequencies of a free ring are compared v/ith exact results 
of Hoppe [l6] and finite difference results of Ball [l9] 
and finite element results of Sabir and Ashwell [34] in 
Table 6,12, An 180^-sector with radial displacement zero 
at the ends was used to find the frequencies of a free- 
ring, from symmetry considerations. 

The results from the present study converge monoto- 
nically from the higher side to exact values. The accuracy 
is considerably higher even for 4 elements 180°-sector 
(13 nodal parameters) results than seven elements 90^^- 
sector (19 nodal parameters^ results of Sabir and Ashwell[343, 
It is interesting to note that this improvement in accuracy 
is achieved despite the fact that the actual element size 
in the present study is almost four times of that used by 
Sabir and Ashwell [34], When the element size is about 
double the size of the element of Sabir and Ashwell [34], 
i,e., for six elements, the results converge to almost 
exact values. The finite difference results [l9] are less 
than the exact values. 

Natural frequencies, of different sectors with 

free—, hinged-, and clamped- end conditions, as shown in 
Figure 6.6 are given in Table 6.13, 
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Table 6.13 


Natural frequencies, of different sectors 

with free-, hinged-, fixed-ends 


Sector 

Angle 

End 

condi- 

tions 

Mode,n 

Natural frequency 

, 10 ^, for 

2 elements 

4 elements 

6 elements 

o 

o 

F„F 

1 jf ^ 3 

0.0 

0.0 

0.0 



4 

8.4107 

8.3914 

8.3910 



5 

24.1399 

23.9376 

23.9267 



6 

49.0616 

47.9738 

47.8124 


S-S 

1 

13.7678 

13.7646 

13.7638 



2 

33,5448 

32.4232 

32.4055 



3 

78.6954 

61.6909 

61.6952 



4 

134,7429 

97.6988 

96.5655 



5 

210.5188 

146.2239 

141.7034 



6 


205.2581 

194.1875 


c-c 

1 

22*6650 

22.6300. 

22.6255 



2 

45*2531 

43,3313 

43.2614 



3 

119*2643 

78.3515 

78.2903 



4 


117.0949 

115.7615 



5 


172.3631 

166.3525 



6 


231.8083 

221.4385 

180° 

F-F 

1,2,3 

0.0 

0.0 

o 

• 

o 



4 

1,8455 

1.8373 

1.8372 



5 

5.3682 

5.3118 

5.3089 



6 

11,4418 

11.1610 

11.1182 


contd. « « 
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Sector 

Angle 

End 

condi- 

tions 

Mode,n 

Natural 

frequency, 


2 elements 

4 elements 

6 elements 


S-S 

1 

2.2679 

2,2669 

2.2667 



2 

7.2115 

6.9280 

6.9238 



3 

18.0850 

13.9824 

13.9828 



4 

32,1272 

23.1352 

22.8484 



5 

50.7828 

35.0758 

33.9399 



6 


50.0236 

47.2303 


c-c 

1 

4.3934 

4.3857 

4.3845 



2 

10.1419 

9.6708 

9.6533 



3 

28.0947 

17.9433 

17.9348 



4 

- 

27.9365 

27.6035 



5 


41.3444 

39.8456 



6 

- 

56.6120 

53.9585 

270° 

F-F 

1,2,3 

0.0 

0.0 

0.0 



4 

0.7699 

0.7601 

0.7600 



5 

2.0374 

2.0021 

2.0006 



6 

4.5613 

4.4294 

4,4082 


S-S 

1 

0.4749 

0.4744 

0.4743 



2 

2.4972 

2.3680 

2.3660 



3 

7.1354 

5.3512 

5.3510 



4 

13.2730 

9.4099 

9.2798 



5 

21.4941 

14.6522 

14.1440 



6 

- 

21.3304 

20.0730 


c-c 

1 

1.4030 

1.3957 

1.3950 



2 

3.7908 

3.5946 

3.5859 



3 

11.4990 

7.0499 

7,0471 



4 


11,5209 

11.3729 



5 


17.3509 

16,6853 



6 

- 

24.2311 

23.0185 


contd , . • 
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Sector 

Angle 

End 

condi- 

tions 

Mode,n 

Natural 

frequency, 

for 

2 elements 

4 elements 

6 elements 

359° 

F-F 

1,2,3 

0.0 

0.0 

0.0 



4 

0.4525 

0.4399 

0,4398 



5 

0.9995 

0.9604 

0.9592 



6 

2.2502 

2.1697 

2.1551 


S-S 

1 

0.0118 

0.0065 

0.0032 



2 

0.9900 

0.9120 

0.9107 



3 

3.4809 

2.4650 

2.4529 




6.8363 

4.7018 

4.5685 



5 

11.4796 

7.6295 

8 .0841 



6 

- 

11.3680 

9.9833 


c-c 

1 

0.5832 

0.5720 

0.5714 



2 

1.7137 

1.6157 

1.6092 



3 

5.8855 

3.4155 

3.4136 



4 


5.8981 

5.8142 



5 

— 

9.1324 

8.7564 



6 

- 

13.0313 

12.3304 


TablG 6,13 shows that convergence is equally efficient 
for all sector angles contrary to the problems encountered and 
mentioned by Sabir and Ashwell [34], namely, the number of 
elements required for large angle sectors is more than that 
for the small angle sectors for the same accuracy. Thus the 
present formulation is simple and more efficient as compared 

to the earlier models. 
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Natural frequencies, for a full ring on different 

number of equi-spaced radial supports. Figure 6.7, are 
presented in Table 6.14. The results are compared with 
the exact values given by Mallik and Mead [24]. The 
accuracy of the exact results [24] is limited by the 
accuracy with which the graph can be read. 

Table 6,14 

Natural frequencies, of a ring on different 
number of equi-spaced radial supports 


No .of 
supports 

Mode, 

n 

Exact 

results 

[24] 

Finite element results, 
elements in each bay 

2 

4 

6 

2 

1 

2.26 

2.2679 

2.2669 

2,2667 


2 

4.38 

4.3933 

4.3857 

4.3845 


3 

6.92 

7.2115 

6.9280 

6.9238 


4 

9.64 

10. 1419 

9.6709 

9.6538 


5 

- 

18,0850 

13.9824 

13.9828 


6 

- 

28.0947 

17.9434 

17.9348 

3 

1 

6.92 

6.9291 

6.9272 

6.9268 


2 

10.44 

10.4291 

10.4102 

10.4086 


3 

10.44 

10.4291 

10.4102 

10.4087 


4 

19.10 

19.8289 

19.0907 

19.0768 


5 

19.10 

19.8289 

19.0907 

19.0768 


6 

23.70 

24.7323 

23.6548 

23.6155 





contd. 

♦ ♦, 
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No. of 
supports 

Mode, 

n 

Exact 

results 

[24] 

Finite element results, 
elements in each bay 

2 ■ 

4 

6 

4 

1 

13.76 

13.7678 

13.7646 

12.3595 


2 

17,90 

17.8948 

17.8724 

17.3260 


3 

17.90 

17.8948 

17.8724 

17.8699 


4 

22.62 

22.6650 

22.6300 

22.6256 


5 

32.28 

33.5448 

32.4233 

33.4557 


6 

37.70 

39.2632 

37.7295 

37.6961 


A good consistent agreement between the finite element 
results and the exact results [24] is observed. The usual 
convergence from the higher side to the exact values is 
obtained. 

The natural frequencies,*^^, of a rotating ring, shown 
in Figure 6.8, for different speeds, ^ , are given in 
Table 6,15. 

The natural frequencies increase with increase in 
rotational speed, 'J5 • The effect of speed is felt more 
on ring with fewer number of supports and on lower modes. 

With increase in number of supports or for higher modes the 
increase in frequency due to rotational speed is marginal. 
This is because as the number of supports increase or it 
vibrates in higher modes the ring becomes more and more 
stiff and, hence, the relative contribution of the rotational 
speed decreases. The inclusion of the rotational speed 
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Table 6,15 

Natural frequencies, of a rotating ring for 
different rotational speeds, li 


No .of 
supp- 
orts 

Rota- 

tional 

speed, 

iX 

Mode, 

n 


Finite element results 

for 

2 ele- 
ments 

4 ele- 
ments 

6 ele- 
ments 

Zero speed 
6 elements 

2 

5.0 

1 

2.8275 

2.8265 

2.8264 

2.2667 



2 

4.7343 

4.7284 

4.7274 

4.3845 



3 

7.4952 

7,2221 

7.2181 

6.9238 



4 

10.3519 

9.8892 

9.8721 

9.6533 



5 

18.2060 

14.1400 

14.1403 

13.9828 



6 

28,1750 

18.0649 

18.0566 

17.9348 


10.0 

1 

3,2932 

3.2920 

3.2919 

Same as 



2 

5.0521 

5.0476 

5.0467 

above 



3 

7.7686 

7,5047 

7.5008 




4 

10.5576 

10.1028 

10.0862 




5 

18.3264 

14.2960 

14.2962 




6 

28.2552 

18.1859 

18.1778 


3 

5.0 

1 

7.1944 

7.1926 

7.1921 

6.9268 



2 

10,6146 

10.5963 

10.5948 

10.4086 



3 

10.6146 

10.5963 

10.5948 

10.4087 



4 

19.9446 

19.2107 

19.1969 

19.0768 



5 

19.9446 

19.2107 

19.1969 

19.0768 



6 

24.8264 

23.7528 

23.7138 

23.6155 


10,0 

1 

7,4503 

7.4484 

7.4481 

Same 



2 

10.7969 

10.7792 

10.7777 

as above 


f 

3 

10.7969 

10.7792 

10,7777 




4 

20.0597 

19.3300 

19.3162 




5 

20.0597 

19,3300 

19.3162 




6 

24.9201 

23.8504 

23.8116 



contd , • . 
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No .of 
supp- 
orts 

Rota- 

tional 

speed, 

-2 

Mode, 

n 


Finite element results for 

2 ele- 
ments 

4 ele- 
ments 

6 ele- 
ments 

Zero speed 

6 elements 

4 

5.0 

1 

13.9198 

13.9166 

12.5281 

12.3595 



2 

18.0141 

17.9919 

17,4489 

17,3260 



3 

18.0141 

17.9919 

17.4489 

17,8699 



4 

22.7605 

22.7258 

22,7214 

22.6256 



5 

33,6158 

32.4966 

33. 5268 

33.4557 



6 

39.3242 

37.7927 

37.7594 

37.6961 


10.0 

1 

14.0702 

14.0670 

12.6946 

Same as 



2 

18.1326 

18.1107 

17.5710 

above 



3 

18 . 1326 

18.1107 

18.1081 




4 

22.8557 

22.8213 

22.8169 




5 

33.6866 

32.5697 

33.5977 




6 

39,3849 

37.8557 

37.8226 



does not pose any additional difficulty and the efficiency of 
the formulation remains uneffected. 

The versatility of the formulation is further illustrated 
by obtaining the natural frequencies, of an I80°~sector 

with three radial supports, placed at unsymraetric positions, 
as shown in Figure 6,9. The results are given in Table 6.16. 
The middle support is placed at varying angles from one end. 
All results are for six-elements in each bay. 
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Table 6.16 

Natural frequencies,"^, of an 180°-sector on 
3-supports with unequal bays 


Mode, 

n 

> 


Angle & 




0° 

rO 

15 

30° 

45° 

60° 

75° 

90° 

1 

2,2667 

4.0121 

5.0194 

6.3908 

8.3085 

11.0282 

13.7638 

2 

6.9238 

9.8697 

12.0082 

14.8978 

18.8429 

22.0896 

17.8705 

3 

13.9828 

18.8553 

22.7949 

28.0826 

33.7344 

27.0225 

32.4055 

4 

,22.8484 

29.6849 

35.6965 

43.5798 

40.3728 

45.7707 

37.6956 


It can be observed that the fundamental frequency increases 
monotonically as the middle support moves from one end towards 
the symmetric position. This gives a significant conclusion 
that the fundamental frequency can be adjusted, by moving the 
middle support position, anywhere between the minimum and 
maximum values for that configuration. The same analysis can 
be extended to rings and sectors on any number of supports, 

6,2.4 Thick Rings 

The natural frequencies, of the free-ring, re = 0.10 
and se = 0.1704, are given in Table 6.17 and compared with the 
analytical values of Kirkhope [36] and experimental results 
of Lincoln and Volterra [IS]* 
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Table 6.17 

Natural frequencies, of a free-ring, re = 0.10 

and se = 0,1704 


Nodal 
dia- 
meters , 
n 

Cl as si— 
cal re- 
sults 
[16] 

Equa- 

tion 

(1.1) 

[36] 

Experi- 

mental 

results 

[18] 

Finite 

2 ele- 
ments 

element 

4 ele- 
ments 

results for 

6 ele- 
ments 

1 

0.0 

0,0 

0.0 

0.0 

0.0 

0.0 

2 

2.6833 

2.5208 

2,5861 

2.5209 

2.5199 

2.5196 

3 

7,5901 

6,6265 

6.8098 

6.8386 

6*5985 

6.5908 

4 

14. 5521 

11.6705 

12.0129 

14.8130 

11.5193 

11.5203 

5 

23.5339 

17.2838 

17.7873 

23.5409 

17.0776 

16.8876 


A good rate of convergence of the finite element results 
can be observed from Table 6,17. As expected, the frequency 
decreases by considering Timoshenko effects in this case also. 
The finite element results are close to the analytical results 
[36], and are on the lower side. Since the finite element 
results are always higher than the exact, the error in Kirkhope’s 
results [36] is more than the finite element results. 

A second comparison of finite element results is made with 
the experimental values of Kuhl [35] and analytical results 
obtained by Kirkhope [36] in Table 6.18. The non-dimensional 
parameters for this ring are re = 0.1378 and se - 0.2015. 
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Table 6,18 

Natural frequencies, of a free ring, re = 0.1378 

and se =*0,2425 


Nodal 

dia- 

meters, 

n 

Classi- 
cal re- 
sults 
[16] 

Equa- 

tion 

(1.1) 

[36] 

Experi- 

mental 

results 

[35] 

Finite 

element results for 

2 ele- 
ments 

4 ele- 
ments 

6 ele- 
ments 

1 

0.0 

0.0 

0.0 

0.0 

0,0 

0.0 

2 

2.6833 

2.4581 

2.4821 

2.4544 

2.4535 

2.4532 

3 

7.5436 

6.2948 

6.1964 

6.4549 

6.2310 

6.2259 

4 

14.5521 

10.8148 

10.4519 

13.5379 

10. 5620 

10.5630 

5 

23.5339 

15.7857 

14.9708 

20.9261 

15.2752 

15.1078 


Consistent with the trend, the frequency decreases further 
with increase in re and se values, the rate of convergence is 
unaffected by this increase. An excellent agreement with the 
experimental values of Kuhl [25] is observed. The analytical 
results of Kirkhope [36] deviate more from experimental values 
compared to the finite element results. 

The ring on radial supports, with Timoshenko effects, 
is analysed and results for the parameters, re - 0.05 and 
se = 0,0866 are given in Table 6.19, The results are compared 

with those given by Rao and oundrarajan [22]. ' 

results, i.e., the results without Timoshenko effects, are the 
finite element results for six elements from Tables 6.12 and 6.1 
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Table 6.19 

Natural frequencies, of a ring on radial supports, 
re = 0.05 and se = 0*0866 


No. of 

Mode, 

Classi- 
cal re- 

Rao and 

Finite 

element re 

suits for 

supp- 

orts 

n 

Sundra- 

rajan 

2 ele- 

4 ele- 

6 elements 

suits , 


6 GlG- 

inonts 

[22] 

ments 

ments 



0 

1 

0.0 

0.0 

0.0 

0.0 

0.0 


2 

2.6833 

2.631 

• 2.6395 

2.6385 

2.6382 


3 

7.5901 

7.292 

7.5699 

7.2979 

7.2918 


4 

14.5583 

13.570 

17.5646 

13.5562 

13.5575 

1 

1 

0.0 



0.0 

0.0 


2 

1.5957 

2.948 

- 

1,5892 

1.5828 


3 

2.4479 

7.740 

- 

2.4127 

2.4126 


4 

5.7745 

14.570 

- 

5.7025 

5.6220 

2 

1 

2.2667 

mm 

2.2397 

2.2388 

2.2386 


2 

4.3845 

3.000 

4. 3266 

4,3181 

4.3169 


3 

6.9238 

7.990 

6,9661 

6.6955 

6.6913 


4 

,9.6533 

15.240 

9.7311 

9.3009 

9.2828 


The results for a free ring match very well v;ith those 

given by Rao and Sundrarajan [22]. All the results are less 

than the classical values, as expected. The results 

rings on radial supports reported by Rao and Sundrarajan [22] 

are inaccurate, because the values are more than the 

, . . „ninn<; The finite element results show the 

ponding classical values, inc 

desired trend and a good rate of convergence. 
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The natural frequencies, of various sectors, 

with the parameters re = 0.05 and se = 0.0866, both ends 
hinged, are given in Table 6.20. The classical results 
are the finite element results for six elements from 
Table 6.13. 


Table 6.20 


Natural 

frequencies, of various 

ends, re = 0.05 and se 

sectors on 

= 0.0866 

hinged 

Sector 

angle 

Mode, 

n 

, Classi" 
cal re- 
sults, 

6 ele- 
ments 

Finite element results for 

2 ele- 
ments 

4 ele- 
ments 

6 element! 

90° 

1 

13.7638 

60.3622 

12,9076 

12.9069 


2 

32.4055 

90.6235 

28.2802 

28.2643 


3 

61.6952 

108.8558 

48.5269 

48. 5306 


4 

96,5655 

190.3344 

69.8358 

69.0849 

180° 

1 

2.2667 

16,5207 

2.2388 

2.2386 


2 

6.9238 

24.1192 

6.6955 

6.6913 


3 

13.9828 

35.1093 

13.0842 

13.0847 


4 

22.8484 

70.7622 

20.8589 

20,6074 

270° 

1 

0.4743 

6.5675 

0.4722 

0.4721 


2 

2.3660 

9.6990 

2.3360 

2,3340 


3 

5.3510 

16.0474 

5.2037 

5.2036 


4 

9.2798 

34.9962 

8.9853 

8.8626 
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The finite element results show a good rate of con- 
vergence and the appropriate trend. The influence of 
Timoshenko effects is felt more on small sectors and higher 
modes, this also confirms to the logic, 

A similar set of results for sector with clamped— 
clamped ends is given in Table 6.21. The classical results 
are again finite elements results for six elements from 
Table 6,13, The non-dimensional parameters are re= 0,05 and 
se = 0.0866. 

These results follow the same pattern as the results 
for hinged-hinged sectors and exactly the same conclusions 
can be drawn. 

The rings on radial supports spinning at constant speed 
are analysed and results are presented in Table 6,22, The 
non-dimensional parameters are re = 0,05 and se = 0,0866. 

The rotational speeds are SJ = 5.0 and 10.0, The six 
element results for a non-rotating ring are also given, 
from Table 6.19, for comparison. 

The trend of the effect of rotational speed on a 
thick ring is same as that on a thin ring. Hence, the 
conclusions from Table 6,22 are same as those from 
Table 6.15, namely, the frequency increases with speed 
and the increase is more for a ring on fewer number of 
radial supports. Here also the convergence rate is not 
affected by the inclusion of rotational speed. 
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Table 6,21 

Natural frequencies, of various fixed-fixed sectors, 

re = 0,05 and se = 0.0866 


Sector 

angle 

Mode, 

n 

Classi- 
cal re- 
sults 

6 ele- 
ments 

Finite 

element results 

for 

2 elements 

4 elements 

6 element! 

90° 

1 

22.6255 

21.0221 

20.9687 

20.9634 


2 

43.2614 

38.5719 

37.2035 

37.1269 


3 

78.2903 

87.6019 

60.7020 

60.5692 


4 

115.7615 

- 

82.3650 

81.4379 

180° 

1 

4.3845 

4. 3266 

4.3181 

4,3169 


2 

9.6533 

9.7311 

9.3009 

9.2828 


3 

17.9348 

25.6310 

16.6959 

16.6826 


4 

27.6035 

- 

25.0164 

24.7265 

270° 

1 

1.3950 

1.3952 

1.3878 

1.3871 


2 

3.5859 

3.7276 

3.5382 

3.5294 


3 

7.0471 

11.0481 

6,8366 

6,8331 


4 

11.3729 

— 

10.9630 

10.8245 


6.2,5 Thin Discs 

The natural frequencies, for a centrally clamped, 

outer-edge free disc, for various values of Ratio = a/b, the 
ratio of inner to outer diameter, are obtained and compared 
with the exact values given by Southwell [4l] wherever possible. 
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Table 6.22 

Natural frequencies, 7 ^^, of rotating rings on radial 
supports, re = 0.05 and se = 0.0866 


No. of 
supports 

Rota- 

tional 

speed, 

-2 

Q 

Mode, 

n 


Finite element results for 

2 ele- 
ments 

4 ele- 
ments 

6 ele- 
ments 

Zero 

speed 

2 

5.0 

1 

2.7924 

2.7915 

2.7914 

2.2386 



2 

4.6623 

4.6556 

4.6546 

4,3169 



3 

7.2402 

6.9797 

6.9758 

6.6913 



4 

9.9326 

9.5109 

9.4933 

9.2828 


10.0 

1 

3.2523 

3.2511 

3.2510 

Same 



2 

^4.9754 

4,9701 

4.9692 

as 



3 

7.5042 

7.2528 

7.2490 

above 



4 

10.1300 

9.7164 

9.6992 


3 

5.0 

1 

6.9552 

6.9535 

6.9531 

6.6965 



2 

10.2040 

10.1829 

10.1812 

10.0022 



3 

10.2040 

10.1829 

10.1812 

10,0022 



4 

18.3661 

17.7129 

17,6993 

17.6887 


10.0 

1 

7.2026 

7.2009 

7.2005 

Same 



2 

10.3793 

10.3587 

10.3571 

as 



3 

10.3793 

10.3587 

10.3571 

above 



4 

18.4720 

17.8228 

17.8093 
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The results are given in Tables 6.23 and 6,24. The value 
of the Poisson’s ratio, v, has been used as 0.3. In case 
of a disc vibrating in n-nodal diameter mode, the natural 
boundary conditions at free end are [60] 


d f r, d^w^ l+2n^ dw . Vn^ dw . (3-v)n^w 
dR ^ ^ dR + ~ dR + ^2 


( 6 . 11 ) 


R 


d^w 

— w 
dR"^ 


+ V 


dw 

dR 


vn^w 

R 


After assembly of equations (5.11) 

( n ) 

end node in the vector ' are 


( 6 . 12 ) 

the terms at the free- 



l+2n^ dw -v 
Rj+r dr'^R=h 


(6.13) 


Bnd = b 

where ND, is the total number of nodal parameters. 

Using the boundary conditions, equations (6,ll) and (6.12), 
the terms of equations (6.13) and (6,14) reduce to 


(3“ '^)n Wj^ 

Bnd- 1 ^dr^R=b 


(6.15) 


where, 



subscript, 



/dw\ 


means evaluated at radius b. 


(6.16) 


These terms must be transferred to the left hand side and 
added at the appropriate places in the stiffness matrix to 
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satisfy the natural boundary conditions. To illustrate 
the influence of these terms on convergence the results 
with improper boundary conditions for Ratio = 0.276 are 
also given in Table 6,23. 

Since the results of Southwell [4l] were obtained 
by evaluating the determinant and then plotting a graph 
to find the root, the decimal accuracy cannot be expected 
i.e,, the accuracy may be limited to the first decimal 
place. A fairly good agreement between the exact and 
finite element results is observed. As can be seen from 
the Table 6,23, the results with improper boundary condi- 
tions are not reliable, because the values converge some- 
times to higher values sometimes to lower values, and the 
trend with increasing nodal diameters is not correct. The 
values increase continuously, instead of decreasing for 
one nodal diameter and then increasing. The rate of 
convergence of results with proper boundary conditions is 
quite fast. 

It can be observed that within the accuracy of 
South§/well ’ s results [4l], the finite element results match 
very 'well. The results converge rapidly for lower as well 
as higher values of Ratio, and for lower and higher nodal 
diameter vibration modes. 
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Table 6,24 

Natural frequencies, of a centrally clamped 

disc vibrating in zero-nodal circle mode for 
various values of ratio 


Ratio 

Nodal 

dia- 

meters 

South- 

well 

[41] 

Finite 

elements results for 

2 elements 

4 elements 

6 elements 

Q.642 

0 

26.21 

26.9426 

26.9174 

26,9162 

0.397 

1 

9.43 

9.4128 

9.3794 

9.3769 

0.603 

1 

22.18 

22.2221 

22.1955 

22.1939 

0.186 

2 

6.55 

6.6246 

6.5644 

6.5514 

0.349 

2 

9.43 

9.4125 

, 9.3538 

9.3473 

0.430 

3 

16.77 

16.4413 

16.3779 

16.3709 

0.590 

3 

26.21 

26.0301 

25.9718 

25.9675 


The results for a centrally clamped disc. Ratio — 0,3, 
are given in Table 6,25, A comparison with the results given 
by Leissa [66] and Mota Soares et al [53] is also made. 

It can be observed that the current results match very well 
with the earlier reported results. A good rate of convergence 
is there. In the current analysis the number of nodal parameters 
for six elements is twleve, which is far less than those reported 
by Mota Soares et al [53j. The accuracy remains of the same 
order. So the use of complicated higher order interpolation 
functions, as in [53], is unnecessary for the problem under 


consideration 
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Tables 6.23 and 6,25 show that the frequency for 
zero-nodal circle modes of disc is less for one-nodal 
diameter compared to the zero— nodal diameter mode of vi- 
bration. For higher nodal circle modes the frequencies 
increase monotonicaily with increase of nodal diameters. 
To study this effect in detail, the natural frequencies, 
tor various values of Ratio, of a centially clamped 
disc vibrating in zero-nodal circle and 0-, 1-, and 2- 
nodal diameter modes are calculated. The results are 
presented in Table 6,26, 


Table 6,26 

Natural frequencies, of a centrally clamped disc 
vibrating in zero-nodal circle modes, for various 

values of Ratio 


Ratio 

Nodal 

diameters 

Finite 

elements results 

for 

2 elements 

4 elements 

6 elements 

.050 

0 

4.2358 

4, 1600 

4.1380 


1 

3.3079 

3.1348 

3.0659 


2 

5.7160 

5.7058 

5.7006 

,100 

0 

4. 5580 

4.4763 

4.4538 


1 

3.8708 

3.7291 

3.6838 


2 

5.9537 

5.9200 

5.9032 


contd 
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Ratio 

Nodal 

diameters 

Finite elements 
results for 


2 elements 

4 elements 

6 elements 

. 150 

0 

4.9776 

4.8973 

4.8789 


1 

4.4749 

4.3530 

4.3197 


2 

6.3194 

6.2622 

6.2383 

.200 

0 

5.5092 

5.4333 

5.4170 


1 

5.1583 

5.0511 

5.0259 


2 

6.8266 

6.7505 

6.7233 

.250 

0 

6.1778 

6.1079 

6.0955 


1 

5.9574 

5.8626 

5.8412 


2 

7., 4999 

7.4107 

7.3850 

.300 

0 

7.0202 

6.9563 

6.9455 


1 

6.9153 

6.8305 

6.8138 


2 

8.3751 

8.2791 

8.2537 

.350 

0 

8.0892 

8.0310 

8.0234 


1 

8.0890 

8.0125 

8.0004 


2 

9.5039 

9.4059 

9.3860 

.355 

0 

8.2113 

8.1536 

8.1467 


1 

8.2211 

8.1460 

8.1310 


2 

9.6334 

9.5351 

9.5129 

.358 

0 

8.2860 

8.2285 

8.2195 


1 

8.3017 

8.2275 

8.2154 


2 

9.7127 

9.6147 

9.5922 


contd 
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Ra'tio Nodal Finite elements results for 

diameters 




2 elements 

4 elements 

6 elements 

.359 

0 

8.3112 

8.2539 

8.2466 


1 

8.3288 

8.2542 

8.2459 


2 

9.7394 

9.6412 

9.6199 

. 360 

0 

8.3364 

8.2793 

8.2705 


1 

8.3561 

8.2817 

8.2675 


2 

9.7662 

9.6679 

9.6477 

.361 

0 

8.3619 

8.3049 

8.2961 


1 

8.3835 

8 . 3092 

8.2981 


2 

9.7932 

9.6948 

9.6736 

.362 

0 

8.3874 

8.3305 

8.3234 


1 

8.4110 

8.3373 

8.3263 


2 

9.8203 

9.7222 

9.7019 

.363 

0 

8.4131 

8.3563 

8,3482 


1 

8.4387 

8.3647 

8.3484 


2 

9.8476 

9.7493 

9.7286 

.364 

0 

8.4389 

8.3822 

8.3746 


1 

8.4664 

8.3929 

8.3815 


2 

9.8749 

9.7768 

9.7545 

.365 

0 

8.4648 

8.4082 

8.4007 


1 

8.4943 

8.4206 

8.4065 


2 

9.9025 

9.8045 

9.7850 


contd 
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Ratio 

Nodal 

diameters 

Finite elements results for 

2 elements 

4 elements 6 elements 

.400 

0 

9.4619 

9.4086 

9.4017 


1 

9.5581 

9.4895 

9.4800 


2 

10.9616 

10.8650 

10.8473 

.450 

0 

11.2526 

11.2040 

11.1969 


1 

11.4387 

11.3776 

11.3790 


2 

12.8607 

12.7678 

12.7492 


From Table 6«26 the change in response of the disc with 
change of the Ratio can be noted very clearly. V^hen the Ratio 
is less than or equal to 0.360, the frequencies for 1-nodal 
diameter and zero-nodal circle are less than those for zero- 
nodal diameter. As the Ratio attains a value more than 0.360, 
the values start increasing monotonically . This indicates that 
for discs having Ratio 5 0.360, the critical frequency is for 

zero— nodal circle, 1-nodal diameter mode instead of the so- 
called fundamental mode of zero-nodal circle and zero-nodal 
diameter. This Ratio of 0.360 is termed as the critical Ratio. 

The effects of rotational speed on the natural frequencies 
of a centrally clamped disc are studied, the results are 
presented in Table 6,27. The results for zero rotational 
speed are the 6-element results from Table 6.23, 


Natural frequencies ^ , of a centrally clamped rotating disc, Ratio = 0,276 
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Rotatio- Nodal Nodal Finite element results for 
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A study of the Table 6.27 reveals some interesting and 
important trends. The same trends are confirmed for various 
values of the Ratio and rotational speeds in Figures 6.10 to 
6.12. As expected the so-called fundamental frequency 
increases monotonically with increase in rotational speed. 

The frequency of all other modes, except the zero-nodal 
circle modes, also increases monotonically, though the 
increase is marginal for higher modes. The critical trend 
is for the frequencies of zero-nodal circles. As the rota- 
tional speed increases the decrease of frequency continues 
even for nodal diameters more than one. The number of 
nodal diameters for which the frequency decreases, increases 
with increase in rotational speed. At higher speeds the 
lowest frequency can be less than the lowest frequency of the 
non-rotating disc, thus, rendering a critical frequency value 
to a higher mode, instead of the so-called fundamental mode. 
Hence, from the design safety point-of-view, the lowest 
frequency must be evaluated. The reason for this decrease 
is the dominating effect of the tangential stress, in 

the 2— D state of stress induced in the disc by the rotation. 
The contribution to stiffness by Ogg is negative, equa- 
tion (5.10). 

A tapered disc has been idealized by constant thickness 
elements, as shown in Figure 6.13. The results for taper 
ratio, T = 0.90 and 0.70 are given in Table 6.28, for a 
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cGntpally clampsd disc. Th© rasults for uniforni disc ar© 
six element results from Table 6.23. 

For a tapered disc the convergence trends are not as 
expected. Table 6,28. The reason for this unexpected trend 
is that, by changing number of elements vie are changing the 
physical system being analysed, as is clear from Figure 6.13. 
Under such conditions the convergence may or may not have 
a definite trend. Hence, for all such cases the element 
matrices must be evaluated by considering the thickness 
variations to have more reliable results, 

6.2.6 Thick Discs 

The numerical results are obtained for a centrally 
clamped, outer-edge free, rotating and non-rotating disc 
including Timoshenko effects, to illustrate the applica- 
bility of the same finite element model for thick discs also. 
The boundary conditions at free-end are satisfied in the 
same fashion as in case of a thin disc, the boundary condi- 
tions are given in Appendix C. 

The natural frequencies, of a centrally clamped 

disc with Ratio = 0.276, and re = 0.02, se = 0,02 are given 
in Table 6.29. The classical results are the six element 

results from Table 6.23. 



Table 6.29 

Natural frequencies, of a centrally clamped disc. Ratio = 0.276 

re = 0*02 and se = 0.02 
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As expected the frequencies decrease with the Timoshenko 
effects, the effect of re and se is more on higher modes. A 
good rate of convergence is observed even with Timoshenko 
effects included. 

The natural frequencies, of a centrally clamped disc 
with Ratio = 0.276 and re = 0.08, se = 0.08 are given in 
Table 6.30. The classical results are the six element 
results from Table 6,23. 


Table 6.30 

Natural frequencies, of a centrally clamped disc, 
Ratio = 0,276, ro = 0,08 and se = 0,08 


Nodal 

dia- 

meters 

Nodal 

circles 

Classi- 
cal re- 
sults , 

6 elements 


Finite elements 
results for 

2 elements 

• 4 elements 

6 elements 

0 

0 

6,5085 

6.3226 

6.2664 

6.2577 


1 

41.6409 

33.2700 

31.4267 

31.2565 


2 

121.1163 

105.5068 

76.8795 

75.4227 


3 

240.9288 

159.7074 

142.6501 

127.2207 

1 

0 

6.3166 

6.1121 

6.0348 

6.0159 


1 

43.8055 

34.4758 

32.7532 

32.5878 


2 

123.6200 

105.7663 

77.7201 

76.3193 


3 

243.7922 

161.4257 

142,9210 

127.8302 

2 

0 

7.7922 

7.3259 

7.2458 

7.2248 


1 

50.6592 

38.1984 

36.8723 

36.7353 


2 

131.7073 

106.7731 

80.4683 

79.2612 


3 

252.5486 

166.2195 

144.0323 

129.9035 


contd,. . 
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Nodal 

dia- 

meters 

Nodal 

circles 

Classi- 
cal re- 
sults, 

6 ele- 
ments 

Finite 

elements results for 


2 elements 

4 elements 

6 elements 


3 

0 

13,5847 

12.0849 

12,0564 

12.0499 



1 

62.8122 

44.5675 

43.8599 

43.7779 



2 

145,7106 

108.9957 

85.5436 

84.6819 



3 

267,5644 

173.3067 

146.4855 

133.8605 



Consistent with the trend the frequencies decrease further 
with increase in values of re and se . Table 6,30 shows a good 
rate of convergence even for higher values of re and se. This 
proves the applicability of the present formulation for the analysis 
of discs including Timoshenko effects. 

To study the influence of rotational speed, ^ , on the natural 

frequencies, a centrally clamped disc with Ratio = 0,276, re = 0,02 

“2 

and se = 0.02 is analysed for rotational speed, « = 5.0 and 10.0. 

The results are presented in Table 6,31, the zero speed results are 
six element results from Table 6,29, 

The trend of the results in Table 6.31 is exactly same as in 
Table 6.27 for thin discs. The only difference is that the values 
in Table 6.31 are lower than the corresponding value in Table 6.27, 
which is as per the expectations. Hence, for thick disc also the 
rotational speed has exactly the same effect as on a thin disc, and 
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Table 6.31 


Natural frequencies, of a centrally clamped 
rotating disc. Ratio = 0.276, re = 0.02 and 

se = 0,02 


Rota- 

tional 

speed, 

w«r 0 

JL 

Nodal 

dia- 

meters 

Nodal 

cir- 

cles 


Finite element results 

for 

2 ele~ 
ments 

4 ele- 
ments 

6 ele- 
ments 

Zero 

speed 

5.0 

0 

0 

10.3607 

10.2303 

10.2032 

6 . 4968 



1 

49.5560 

46,1021 

45.7292 

40.7332 



2 

185.9770 

125.0171 

120.5761 

115.4480 



3 

620,5471 

263.1034 

228.8881 

223.8489 


1 

0 

9.4644 

9.2870 

9.2456 

6.3013 



1 

51.0212 

47.6286 

47.2532 

42.8077 



2 

187.2666 

126.8943 

122.5885 

117.7245 



3 

621.2715 

264,7450 

230.9235 

226.0236 


2 

0 

8.0768 

7.8052 

7.7356 

7.7684 



1 

55.7588 

52.6977 

52.3469 

49.3448 



2 

191.2148 

132.7407 

128.8828 

124.7605 



3 

623.4758 

269.7568 

237.1433 

232.6288 


3 

0 

10.8020 

10.6400 

10.5905 

13.4788 



1 

64.5425 

62.2994 

62.0424 

60.8269 



2 

198.0335 

143.0786 

140.0364 

136.9882 



3 

627.1676 

278.3696 

247.8678 

243.9486 


contd 
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Rota- 

tional 

speed, 

* Z.I 

Q 

Nodal 

dia- 

meters 

Nodal 

circles 


Finite 

element results 

for 

2 ele- 
ments 

4 ele- 
ments 

6 ele- 
ments 

Zero 

speed 

10.0 

0 

0 

13.0641 

12.8717 

12.8263 




1 

53.7715 

50.6083 

50.2235 




2 

188.9033 

129.6675 

125.4605 




3 

621.7218 

267.1426 

233.8097 



1 

0 

11.6969 

11.4333 

11.3595 

Same as 



1 

54,8849 

51.7122 

51.3170 

above 



2 

190.1133 

131,3409 

127.2344 




3 

622.4418 

268,6986 

235.7195 



2 

0 

8.0816 

7.5599 

7.4255 




1 

58.6385 

55.6045 

55.2125 




2 

193.8299 

136.6070 

132.8537 




3 

624.6108 

273.4569 

241.5722 



3 

0 

6.6061 

6.0114 

5.8397 




1 

66.0247 

63,5953 

63.2814 




2 

200.2841 

146.0698 

143.0113 




3 

628.2533 

281.6684 

251.7189 



same conclusions can be drawn. A good rate of convergence can be 


observed from Table 6.31 
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6 . 3 SUMMARY 

The finite element formulations for the beams, the 
rings, the sectors, and the discs including Timoshenko 
effectsh-avo been proposed in a systematic and logical 
fashion# The models proposed are the simplest requiring 
minimum number of nodal parameters per element. The 
formulations proposed are capable of analysing the effects 
of rotational speed or prestressing. The efficiency and 
accuracy of the formulations is illustrated. 

The effects of ratio of inner diameter to outer dia- 
meter and rotational speed on vibration modes of a disc 
are studied very critically. Two important conclusions 
have been drawn in case of the discs. 



CHAPTER VII 


CONCLUSIONS 

The following conclusions axe drawn from "the present 

study : 

1* Systematically and logically developed finite element 
formulations for the structural components, beams, 
rings, sectors and discs, including Timoshenko effects, 
are simple and efficient compared to all the earlier 
proposed formulations. The efficiency of the formula- 
tions is illustrated by many numerical examples. The 
accuracy is proved by comparison of the results with 
exact/experimental/approximate results wherever 
possible. 

2. The inclusion of the Timoshenko effects, rotational 
speed, pre-stress, etc,, in the analysis of the stru- 
ctural components does not necessitate the change of 
interpolating function i.e,, the interpolation poly- 
nomial remains same irrespective of the fact that 
whether these effects are included in the analysis 

or not* 

3. In case of an 180°-sector on two supports, if a third 
support is added at varying angles from one support, the 
fundamental freguency increases monotonically with 
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increase in the angle till the symmetric position is 
attained. A similar trend can be expected for other 
configurations. 

4* In a non-rotating centrally clamped disc with inner to 
outer diameter ratio less than or equal to 0.360, the 
minimum frequency occurs for zero-nodal circle, one- 
nodal diameter mode instead of the zero-nodal circle, 
zero— nodal diameter mode. As the inner to outer dia- 
meter ratio becomes more than 0.360 the minimum fre- 
quency is for the zero-nodal circle, zero-nodal dia- 
meter mode, 

5. Depending upon the rotational speed, the lowest frequency 
in a centrally clamped disc may occur for any mode with 
zero-nodal circle and any nodal diameter, irrespective 
of the value of the inner to outer diameter ratio. 
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APPENDIX A 


BOUNDARY CONDITIONS AInID ELEMENT MATRICES FOR 

BEAMS 

A.l BOUrOARY CONDITIONS 

The boundary conditions for the Timoshenko beams are 
[3,55] : 

i) Hinged end 

w = o (A.l) 

■ijj' ss 0 (at extreme boundaries only) (A. 2) 

ii) Clamped end 

w = o (A. 3) 

^=0 (A. 4) 

iii) Free end 

' = 0 ( A. 5) 

i w’- '1' = 0 (a. 6) 

where L is the length of the beam. 

A. 2 ELEMENT MATRICES 

The various constitutive matrices for the elemental 
stiffness and mass matrices are given below. 
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A. 2.1 Approach A 


[A^] = I {N” } LN"J dx = -i- 
0 h" 


12 

symmetric 

6h 

4h^ 


-12 

-6h 

12 

6h 

2 

2h^ 

-6h 4h^ 


[A,] = /{N' ILN’J dx =1 


[Aj = /{N } LNJ dx 


h/lO 2h"^/l5 symmetric 
-6/ 5 »h/lO 6/5 

h/lO -h^/SO -h/lO 2h^, 


156 

2 

22h 4h symmetric 
54 13 h 156 

-13h -3h^ -22h 4h^ 


( El - (PI+EI 


{Y,} 


33je) P 

^ 0 / 


{y^ } * 




3 X Ij 
0 
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where subscripts j and k means the evaluation at jth and kth 
nodes respectively. 


A. 2. 2 Approach B 


The linear variation of the cross-sectional properties 
can be expressed as 




and 


.(e) 




/• N 

I . 

•k 




The constitutive matrices for this variation are : 
h 




“ aG(A.+A^) 

i aAG (N'} lN’J dx = — ^ 


o 


L 


1 ”1 

1 1 


[Ko3 


^ aAG ( N' } L^^Jdx = ^ 


2Aj4'Aj^ Aj4*2Aj^ 

-(.2A.+Aj^) -(A.+2Aj^) 


[K 3 / (El ( N» Un’ J + aAG (N > LNJ)dx = 


— 


h 




- 


1 -1 

, aGh 

3Aj+A^ 

A^+A|j 

-1 1 

+ 


Aj+3\_ 


i P In 

[Mj = i PA{N )LNJ dx 


0 


3Aj+Ak Aj+Aj^ 
Aj+Aj, Aj+3Aj^ 
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[Mo] 


/ PI (N 1 LNJ dx = -jI 


{B2> 


(ne) 


3IH-I, 




I..3I, 


■( aAG( 
( aAG( 


3w(®> 

,(e) 

3x 




3x 


(El ir~) k 


such that the elemental vector of nodal parameters is 


{ w 1 



[Ko] matrix using l~point Gauss-Legendre integration formula is 


Ed.+Iu) 

[K,] = — ifT^ 


1 ”1 

•1 1 


(xG( A^+Aj^)h 

+ f 


1 1 
1 1 



APPENDIX B 


BOUI^ARY CONDITIONS AND ELEMENT MATRICES FOR RINGS 

B . 1 BOUNDARY CONDITIONS 

The boundary conditions for rings are [57] : 

i) Pinned end 


w = 0 


9w 

30 


= 0 


3^w 


0 


ii) Clamped end 


w = 0 


3w 


= 0 


30 ^ 


0 


(B.l) 

(B.2) 

(B.3) 

(B.4) 

(B.5) 

(B.6) 


iii) Free end 


3w 

30 





(B.7) 
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^ 9'^w 
30^ 38"^ 


0 


3w 

36 


3w 

90 


3 © 


= 0 


(B.8) 


(B.9) 


B . 2 ELEMENT MATRICES 

The various constitutive matrices for the elemental 
stiffness and mass matrices are s 




[A^] = / {N»’ ' } LW’Jde = I 


r240 

I20p 


symmetric 


64f 


20p^ I2p^ 

1-240 -I20p -20p^ 240 
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120p 56p' 


i -20p' 
L 

600 


-sp' 
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o 4 


64p' 
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J 
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^ 35B 
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300p 192P' 
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2 
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4 

600 


4p^ -coop I92p^ 
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L. 


J 
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{^2}= 


0 
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0 
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3 0^ 


+ 2 - 
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APPENDIX C 


BOUNDARY CONDITIONS AND ELEfAENT MATRICES FOR DISCS 


C.l BOUNDARY CONDITIONS 

The boundary conditions for a thick disc, vibrating 
in n-nodal diameter mode, are [59] t 


i) Hinged end 

w = 0 

vn j^-0 

ii) Clamped epd 

w = 0 



(C.l) 

(C,2) 

(C.3) 

(C.4) 


(C.5) 


(C.6) 
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C.2 ELEMENT MATRICES 


The stresses induced due to rotation are given 

°RR = ^ 

""ee = =1 + ^ - =4®^ 


where 


= 15 + = 3 '=' 


C„ = i P 

2 8 a^+d+vjb^ 


C3 = ^ ( 3+ v) p 


= I ( 1+3 V ) p 


by [58] 
(C.7) 

(C.8) 

(C.9) 

(C.IO) 

(C.ll) 

(C.12) 


The various constitutive matrices for the elemental 
stiffness and mass matrices are given below. 


[A,] 


h 

/ (R.+r) {N" } LN "jdr = 

r\ J 



12 

6h 

0 

4h symmetric 

1 

6 

2h 

h^ 

symmetric 

-12 

-6h 12 


-6 

-2h 

6 

6h 

2h^ -6h 4h^ _ 

i 

4h 

h^ 

-Ah 
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[A 2 ] = J {N«}LN'jdr = 



36(F2-2F3+F^ 

-6h( F.-5Fo+ 
7 F 3 - 3 F 4 ) ^ 


)• 

h^{ F q- 8 F^+ 

- 22 F 2 - 24 F 3 + 9 F 4 ) 


symmetric 


-36(F2-2F3+F4) 6h(F^- 5 F 2 + 36( F 2 - 2 F 3 +F 4 ) | 

7 F 3 -. 3 F 4 ) I 

6h(2F2-5F3+ ~h^( 2 F^- 11 F 2 + -6h(2F2-5F3+ h^( 4 F 2 --I 2 F 3 +' 

3F4) ■ I8F3-9F4) 3F4) 9F4) 


where F^ = ln(l+h/R.) (C.13) 

1 

Fi = j ^ 6 (C.14) 

h 

["^ 3 ] = / -, ^ ‘ ' , ' 3 tm tNJdr = 

Q (Rj+r) 


Z 0 - 6 Z 2 + 4 Z 3 + 

9 Z 4 -. 12 Z 5 + 4 Z^ 

n symmetric 

h(Z^- 2 Z 2 - 223 + h"^(Z2-4Z3+ 

8Z4-7Z5-h2Z^) 6Z4-4Z5+Z^) 


3 Z 2 - 2 Z 3 - 9 Z 4 + h(3Z3-8Z4+ 

12Z5-4Z^ ^Z^-ZZg^) 

-h ( Z2"'^3~ 8Z4+ -h^( Z3- 3Z4+ 

SZ^-ZZ^) 3Z5-Z^) 


9Z4-12Z5+4Z^ 

-h(3Z4- h^(Z4-2Z3+Z^) 

SZ^+ZZ^) 


where 


h(R.+Rj^) 
— -i-- — 

zti. 


(C.15) 
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(R.+Rj,) 


h ln(Rk/Rj) - 


R4 + R.. 


2h R. 


(C.16) 


(C.17) 


1 3R. 3R. R.(R.+R^) 

i: _1 T^fT>/n\ . _ .1 .1 K-' 


Z3 = — ln(Rj^/R. ) + - 


h h"' ^ J h Rj^ 2h^Rj 

1 , . 4R. 6 R^ 4r2 R?{R.+R^) 

jJ: ^ n I r> \ 1 .1 1 / rv / 7-1 \ 1 1 1 ic 


Z s: R*+R, + -.^Xn(R|./R- ) -f — 

^ 2h'^ J h'^ ^ h\ 2h^Rj 


(C.18) 


R^-R^ 5R. 


lOR? lOR? 


xvjn. 

^ (Rj+Rk)+ ln(Rj^/Rj) + 


5R? f^i(Ri+Rk) 


h^R, 2h'^R? 

K, o *>• 


R^-R^ 6R.(rJ-R^) 15R^ 


15R- 6R^ RT(R 

— ^ ln(R,/R. ) ~ H— * + ■■‘^■s--5 
h'^ K j 2h^R^ 

h 

[A4] = / (J^^(Rj+r) {N» }tN»Jdr 


—4. 4 (R +D ) 4. 

2h^ ^ 

6Rf Rf(R,+R,) 


(C.19) 


(C.20) 


(C.21) 


Co(2 


(C^-R^C3)[A^]-C2[A2]- 


^ symmetric 
h/10 hV30 

-3/5 -h/10 3/5 

0 -h^/60 0 h^/lO 


12/35 


symmetric 


h/14 2hVl05 
+3R-h 10 

J -12/35 -h/14 -§1 

-h/35 -h^/70 ^ 
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3/14 




h/20 

o , 

symmetric 

llhV840 



-3/14 

-h/20 

3/14 


-h/28 

-llh^/840 h/28 

13hVl68 


0 

[Aj] = ^ {N }tNJclr = C [Ag] + CjEA,] - CJA^] 

^ o R.+r ^ ^ ^ ^ 

h 

[Ag] = / (Rj+r) {N ) [N'Jdr = 

o 



6/ 5 


3/5 

^ symmetric 


^ symmetric 

h/lO 2hVl5 

1 

h/lO hV30 

-6/5 -h/10 6/5 

T 

-3/5 -h/lO 3/5 

h/lO -h^/30 -h/lO 2h^/l^ 


0 -h^/60 0 h^lO 


C A-^] 



h 

= / (R.+r) {N }LNJdr = 

o ^ 


156 


3/35 

^ symmetric 

22h 4h"^ 


^ symmetric 

h/60 hV280 

54 13h 156 

+ 

9/140 h/60 2/7 

-13h -hV30 -22h 4h^: 


-h/lO -h^/280 -h/28 h^/l68 
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[ Aq] 


- / LNJdr = 

o j 




h 


9F4-12F5+4F^ 


symmetric 


h ( F j^- 2F2-2F3+ h^{ F2-4F3+ 

8F4-7F5+2F5) 6F4-4F5+FJ) 


h( 3F3-8F4+ 9F4-12F3+4Fg 

’^5-2%) 


Temp 


3F2»2F3~9F44- 
12F5-4F^ 

-h(F2-F3-3F4+ -h^( F3-3F4+ -h( 3F4-5F3+2F^) h^( F4-2F5+I 

_SF5-2F^) SF^-F^) F^)^ 

D -L (R — l+2n^) D i - 


3r 


3r 


(PI D R • - TR 0 ^ 


(e) 




rr 3 r 


{Yi } 


Temp 

< ^ > 
Temp Ij 

0 


{Y^}: 


0 




-D(R 

3r j 


0 


^ (e) 

D(R^ )], 



